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Abstract 

Off-diagonal vacuum and nonvacuum configurations in Einstein gravity can mimic physical effects of 
modified gravitational theories of f(R , T, R llv T tw ) type. To prove this statement, exact and approximate 
solutions are constructed in the review, which encode certain models of covariant Horava type gravity with 
dynamical Lorentz symmetry breaking. The corresponding FLRW cosmological dynamics with possible 
nonholonomic deformations and the reconstruction procedure of certain actions closely related with the 
standard ACDM universe are studied. Off-diagonal generalizations of de Sitter universes are constructed 
which are generated through nonlinear gravitational polarization of fundamental physical constants and 
which model interactions with non-constant exotic fluids and effective matter. The problem of possible 
matter instability for such off-diagonal deformations in (modified) gravity theories is briefly discussed. 

Keywords: Exact solutions of modified gravities, accelerated expansion of the universe, cosmological 
reconstruction procedures, /(i?)—gravities and generalizations, ACDM cosmology. 

PACS numbers: 98.80.-k, 04.50.Kd, 95.36.+X 


1 Introduction 

There are several motivations for the study of modified gravity theories. At short scales it seems clear 
that General Relativity (GR) needs to be modified in order to take into consideration quantum effects, which 
leads to second- and higher-order terms in the curvature R, and these same modifications might be also very 
useful to solve problems at large scales, as the acceleration of the universe expansion and other, as the nature 
itself of dark matter (DM) and dark energy (DE). Not the least, it is an attempt to formulate a self-consistent 
(in some particular way, at least) theory of quantum gravity. Moreover, an increasing amount of more and 
more accurate and constraining observational data will help to discriminate among the different, alternative 
modifications of the gravity theory in the search for a better description of our universe. 
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Among the different classes of extensions of general relativity some of the most popular are f(R), f(R, T ), 
and /(R, T, F )—which we will here generically call /(R,...)-modified theories, being R the Ricci scalar and 
T the metric torsion. In these approaches, the standard Lagrangian for GR, namely as C = R, on a 
pseudo-Riemannian manifold, V —where R is the Ricci scalar curvature for the Levi-Civita connection, V— 
is modified by the addition of a functional, f(R ,...), of the Ricci scalar only, in the first case, of R and the 
torsion tensor, Tp , the energy-momentum tensor for matter, Tg 7 , and/or its trace T = T“ (in the second), 
and of a generalized Ricci scalar R, and a Finsler generating function, F , in the third case (such values may 
be defined on the tangent bundle TV), etc. Classes of modified theories of these kinds can be successfully 
constructed, and also the corresponding reconstruction procedures, able to mimic the ACDM model including 
the dark energy epochs and the transitions between the different main stages of the universe evolution, thus 
providing a unified description of the entire cosmological history. For reviews of some of the most important 
results along this line, see Refs. [Tj, [2j |3j SJ ,5] [6, 7, |8j 9j [TU, Ft, 72; [T3], 77] . 

Several of these model constructions of modified theories are actually related with different forms of the so- 
called covariant Horava gravity associated with a dynamical breaking of Lorentz invariance dang naiig, 
and with further developments, as well, including in particular generic off-diagonal solutions, Lagrange- 
Hamilton-Finsler like generalizations, deformation quantization, A-brane models, and gauge like gravity 
mmmm- In some simplified approaches, theories of this kind can be constructed in a power-counting 
renormalizable form or as nonholonomic brane configurations which correspond to power-law versions of 
actions of type f{R, T, R.^ V T ,W ) [2711251 .26]. In general, the spacetime geometries can be of Finsler type, with 
commutative and/or noncommutative parameters, and off-diagonal metrics for wapred/trapped solutions 
[271 M 1291 EB EH). This also includes effects as Lorentz violations, nonlinear dispersion relations and 
locally anisotropic re-scaling, and effective polarizations of constants, which provide a deeper understanding 
of the possible connections between these f(R ,...) modified theories, and Horava-Lifshitz and Finsler theories 
(see also Refs. [M1E31E1E3IM1E3ESJE5I, in this respect). 

Field equations for gravitational and matter field interactions in GR and various modified theories of the 
types described above usually consist of very sophisticated systems of nonlinear partial differential equations 
(PDEs). No wonder they request advanced numeric, analytic and geometric techniques for constructing 
exact and approximate solutions. The most important physical solutions, for black hole configurations, 
observable cosmological scenarios, etc., have been therefore constructed with the simplifying diagonalizable 
ansatz for the corresponding metric (obtained by appropriate coordinate transformations and frame rotations) 
with Killing symmetries. After a series of assumptions of “high symmetry” of the relevant interactions (for 
spherical, cylindrical or torus ansatze, with a possible additional Lie group interior symmetry), the systems 
of resulting nonlinear PDEs are usually transformed into much more simplified systems of nonlinear ordinary 
differential equations (ODEs), what is a great advantage, indeed. In such cases, some classes of exact 
solutions can be obtained in explicit form (see the monographs maim for reviews of some results in GR). 
Even then, it actually took more than half a Century to understand the fundamental physical implications 
of these solutions; for instance, of the Schwarzschild and Kerr black hole metrics, and to finally elaborate the 
Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmological scenario. At present, it is already possible to 
construct more general classes of exact solutions in (modified) gravity theories depending generically on two 
or three variables in a “less symmetric” form for four-dimensional (4-d), and extra-dimensional models using 
advanced geometric, analytic and numerical methods. We can provide a plausible physical interpretation 
when such solutions are defined by certain symmetry transforms or small deformations of some well-known 
classes of solutions. However, to derive a new class of exact solutions of a gravity theory is by itself not of 
main interest for physicists, neither to applied mathematicians, unless such construction does result in new 
and interesting classical or quantum-physical effects, or does provide convincing explanations to cosmological 
data coming from the newest surveys. 

In a series of works EZlEHlEHlEniEH, the so-called anholonomic frame deformation method (AFDM) 
for the construction of exact solutions in gravity has been developed. It provides a general geometric tech¬ 
nique which allows to integrate PDEs for gravitational and matter fields interactions, for generic off-diagonal 
metrics with generalized connections, or for the torsionless Levi-Civita one. Such solutions may depend on 
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all spacetime coordinates via various classes of generating and integration functions, commutative and non- 
commutative parameters, and so on. In particular, the corresponding metrics and connections can exhibit 
anisotropic ellipsoidal or toroidal symmetries, or encode certain generalized solitonic hierarchies, coming from 
an effective non-holonomic (with non-integrable constraints) dynamics with nontrivial topological configura¬ 
tions. It is furthermore possible to analyze the physical implications of geometric constructions of this kind, 
provided they describe certain “small parameter” deformations related to well defined black hole objects in 
cosmological models, or particle physics interactions with broken symmetries. 

In brief, the AFDM is based on a quite surprising decoupling property of the vacuum and on certain 
classes of non-vacuum fundamental field equations in GR and modified gravities. The main idea is to work 
with an “auxiliary” connection when physically important systems of nonlinear PDE decouple in certain 
classes of nonholonomic frames. This allows to integrate systems of this kind of very general form. Usually, 
the auxiliary connection which is needed involves nontrivial nonholonomically induced torsion, which in GR 
is determined by certain generic off-diagonal terms of the metric and corresponding classes of nonholonomic 
(equivalently, anholonomic, i.e. non-integrable) constraints on gravitational and matter field dynamics. A 
nonholonomically induced torsion is different from that in the Einstein-Cartan or string gravity theory, where 
torsion fields are subject to additional algebraic or dynamical field equations. Having constructed certain 
general integral varieties of solutions, we can almost generically consider certain classes of constraints where 
the auxiliary connection transforms into the Levi-Civita one. Here we note that it is important to impose 
such zero-torsion constraints after certain classes of generalized solutions are found in general form, but 
not before applying the AFDM. Generic off-diagonal solutions for nonlinear systems can be restricted to 
torsionless configurations, provided certain nontrivial solutions have been already found. If non-holonomic 
constraints are imposed from the very beginning (for instance, spherical symmetries and a simplified diagonal 
ansatz for the metric), then one excludes from the analysis more general classes of nonlinear interactions. 

The crucial importance of generic off-diagonal solutions in GR and modified theories is determined 
by a series of geometric and analytic properties of the associated systems of nonlinear PDEs which are 
used for elaborating the cosmological models and for performing the quantization of the associated gravity 
theories. In this respect we should emphasize three key issues: 1) A number of physical effects and observed 
cosmological data can already be explained in GR or in the context of modified gravity theories (for instance, 
within f(R) gravity or with nontrivial massive terms, see pl2l 33] ) or can alternatively be modelled with 
the help of off-diagonal interactions in different types of modified theories. 2) Generic off-diagonal solutions 
encode configurations with a nontrivial parametric vacuum and effective matter field interactions, also with 
gravitational polarizations of the interaction and a cosmological constant, and for nontrivial generating 
and integration functions. For certain well-defined conditions, such models describe broken fundamental 
symmetries—as, for instance, locally anisotropic interactions, terms deploying violation of local Lorentz 
symmetry, and warped and trapping configurations—which points towards new methods of quantization and 
results in alternative scenarios of accelerating and anisotropic cosmological theories, with different solutions 
for the dark energy and dark matter physical problems. 3) Considering off-diagonal configurations one can 
model classical and quantum f(R) modified gravities and the like, and also Horava-Lifshitz and Finsler like 
theories in a suggestive, unified geometric way [121 [201 [2D [221 EZl [2SI [291 SOI El]. 

The aim of this paper is to apply the anholonomic frame deformation method for the constructing of 
exact off-diagonal solutions corresponding to cosmological models of modified gravity of the general form 
f(R, T, R^yT^ w ), and to study the conditions under which such configurations can be alternatively modeled as 
effective Einstein spaces with nontrivial off-diagonal parametric vacuum and non-vacuum configurations. The 
FLRW cosmological dynamics and a reconstruction procedure of the ACDM universe will be investigated. To 
be noted is that we will not work with exotic anisotropic fluid configurations as in nsuisi nzi ns eh Ha E6| , 
but rather with off-diagonal deformations of de Sitter solutions, as in d2U22U2U E2] - The problem of matter 
instabilities in modified and deformed GR theories will be analyzed and solutions will be obtained for certain 
classes of non-holonomic configurations. 
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2 Off-diagonal interactions in modified gravity and cosmology 


In this section, we formulate a geometric approach to f(R , T, R^ IV T^ U ) gravity and summarize the anholo- 
nornic frame deformation method 123 ESI 1221 EDI EH- The geometric constructions will be adapted to non- 
holonomic distributions with associated nonlinear connection (N-connection) structure^ The N-connections 
formalism will be used for constructing certain classes of N-adapted frames with respect to which the grav¬ 
itational and matter field equations decouple in very general forms (see Sect. [3] below). This is possible for 
metric compatible linear connections with non-holonomically induced torsions completely defined by metric 
tensors. Imposing additional constraints, we can generalize zero torsion configurations for the Levi-Civita 
connection, V. 

2.1 Modeling dark energy with off-diagonal metrics 

In order to motivate our approach, we discuss simple FLRW cosmology and dark energy and dark matter 
models which are extended for generic off-diagonal solutions in GR and modifications. Working in a spatially 
flat spacetime with diagonal quadratic form 

ds 2 = g a (t)(du a ) 2 = a 2 (t)[[dx 1 ) 2 + (dx 2 ) 2 + ( dy 3 ) 2 ] — dt 2 , ( 1 ) 

for local coordinates u a = (x*,y 3 ,r / 4 = t), when i = 1,2, the FLRW equations are 

-^H 2 = p and jf + 3 H (p + p) = 0 , 

where p and p are, respectively, the total energy and pressure of a perfect fluid (pressureless or just radiation), 
H := a°/d for a° := da/dt = = dtd, and k 2 is related to the gravitational (Newton) constant!! To explain 

the observational data of an accelerating universe, and dark energy and matter, various models have been 
studied (see reviews and references in ID El El SI El El El El El HQ EUniEl EH), with effective or exotic 
matter with an equation of state (EoS) of phantom kind, p = zup. with w < —1. The simplest model of 
phantom DE is given by 

3 

—Hde = pde and p'de + 3Hde(1 + vj)p£>e = 0 , 
which for w < — 1 admits an exact solution 

2 

dd-DE w, • (^) 

3(1 + w)(t s — t) 

This solution has a finite-time future singularity (Big Rip) at t = t s . 

Some models have been considered where the Hubble function H(t) is determined by a phantom DE 
coupled with DM, via a coupling constant, Q, which results in the conservation law 

Pde + 3H(1 + xu)pde = —Qpde, Pe>m + 3 Hpdm = Qpdm- 
The solutions of these equations can be expressed as 

Pde = °Pde e - 3 ( 1 + ro ' ) e -< ’^ and pdmo 3 = Q°pde J dt' e -3a7 e - ^ , 

1 It should be noted that in generalized Finsler like theories, the N-connection structure is given by a set of three fundamental 
geometric objects which, for certain models, is completely defined by the so-called Lagrange/Finsler generating function. We 
do not study in this work Finsler like modifications of GR. 

2 We use a system of notations different from that in standard cosmology, as this will be convenient for constructing cosmo¬ 

logical models with generic off-diagonal metrics, and also in order to follow the conventions in our previous works. 
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respectively, where °Pde is an integration constant and the EoS is taken to be p = wp^E- These functions 
can be used for the second FLRW equation, 


-4(2#* +3 H 2 )=p. 

We have 

H = -Q/ 3(l + ro) (3) 

for the exact solution of this equation, which corresponds to the evolution for de Sitter space, a(t) = 
aoe - ® f//3 ( 1+ro ), where oo is determined from a 3 l 1+ro ) = — ^r(l + vu) 2 w °pde • The value of H in ([3]) is 
positive for w < —1, what does not mean that the Big Rip singularity in (J2]) can be avoided, but just shows 
that the coupling of the phantom DE and DM gives a possibility that the universe could evolve as a de Sitter 
phase. More than that, the first FLRW equation, 



Pde + Pdm, 


(4) 


imposes the relation pdm = (1 + w)pde- Considering a de Sitter solution as an attractor, with w ~ —4/3, 
we obtain —(1 + w) ~ 1/3, which is almost independent from the initial condition, i.e., it solves for free the 
so-called coincidence problem!! 

Since the DE-DM coupling does not always remove the singularity and there is no such fluid with 
constant EoS parameter, models were considered which are proportional to a power of the scalar curvature, 
for instance, pfluid, oc R l+e , for e > 0. In that case the total EoS parameter is greater than —1 and a Big 
Rip does not occur for large curvature. Two variants of theories have been exploited where this kind of 
inhomogeneous effective fluid matter is realized, by a conformal anomaly and other quantum effects or by 
some modified model of gravity, for instance, when the gravitational Lagrange density R —)■ f(R) = R + R x . 
In the case 1 < x < 2, we have that solutions with 


Ht ~ 


(x — l)(2xr 
x- 2 


1 ) 


and w e ff ~ — 1 — 2iL°/3iL 2 > —1 


do not result in a Big Rip or any other kind of future singularity. Similar classical and quantum arguments 
were considered as motivations to study f{R) modified gravity theories [T| [2, 3j 01 [5, 5j El SI 13 EDI El, HI 

Ha eh ns Ea ezi esi eh esi ehi hh ssj. 

In a series of works EZI C2S1G3 EDI Em E3 GDI EE EH ® S3 03 07] , various classes of off-diagonal solutions 
were studied which can be constructed by geometric methods in modified gravity theories. We proved that, 
for instance, certain important effects and cosmological models related to f{R) modified theories and the 
like can alternatively be explained by nonlinear off-diagonal gravitational and matter held interactions with 
respect to nonholonomic frames. Let us briefly recall the main ideas supporting such an approach. Off- 
diagonal ansatze for metrics (see, for instance, (1591) 1. 

gi + u> 2 {w 2 h^ + riihi) Ld 2 {w\W 2 h^ + 7 x 1712 / 14 ) oj 2 w\h 3 w 2 ni/i 4 
_ oj 2 (wiW2h3 + ninzhi) 52 + w 2 (w 2 2 ^3 + n^hi) w 2 w 2 ^3 u 2 ri2hi /r,, 

UJ 2 W\h3 W 2 W2^3 W 2 /l3 0 ’ ^ ' 

W 2 7li/l4 W 2 ri2^4 0 U> 2 hi 

where the coefficients are parameterized by functions of type g± = §2 ~ eand rii(x k ) (we can fix 
certain constants for corresponding classes of generating, <h(x fc ,t), and integration functions), h a [&(x k ,t)] ~ 
h a (t), [for a = 3,4], Wi[&(x k , f)] ~ Wi(t ) and uj(x k ,t ) ~ u>(t), were found to generate exact (in general, 
nonhomogeneous) cosmological solutions in modified gravity theories. Such generic off-diagonal metric^ can 
be represented in the form 

ds 2 = ^(^[(e 1 ) 2 + (e 2 ) 2 ] + a 2 (f)/r 3 (t)(e 3 ) 2 + (e 4 ) 2 , (6) 

3 If DE does not couple with DM, we have pdm ~ a -3 and pde ~ a -3 ^ 1+ro \ which do not satisfy the observed 1/3 ratio of 
DE and DM and does result in a coincidence problem. 

4 Which cannot be diagonalized by coordinate transformations. 
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with respect to so-called N-adapted frames (Eqs. (11211 and (1131) are used, in general) 

e 3 = dy 3 + riidx l ,e‘ l = dt + Wi(t)dx l . 

For certain well-defined conditions (see Sect. ED, we can consider off-diagonal deformations g a (t) — >• 
9aj3(x k ,t ) r\_/ g a p(t ) defining new classes of cosmological models which mimic contributions from f(R) modi¬ 
fied gravity encoded into the data for uj(t),Wi(t), etc. The corresponding formulas are nonlinear functionals 
relating generating functions to the (effective) matter sources. Such off-diagonal configurations are equiv¬ 
alently modeled as solutions of some effective field equations R°g = A 8 a p. This way, various classes of 
cosmological solutions of modified gravities can be alternatively modeled by metrics of type ©, when the 
scaling factor a(t) is nonlinearly determined by the coefficients Wi(t) and h a {t) via a generating function 
<f>(f) and an effective source T (t). We can model ACDM cosmology and analogously DE and DM effects 
with p£>E + Pdm encoded into <h(t) and T (t), but with respect to the adapted frames e a (t). Solutions with 
off-diagonal metrics may be interpreted in accordance with observational data if a(t) is chosen to determine, 
for instance, an effective H(t) (0 with cosmological evolution from a spacetime background encoding f(R)- 
modifications. To prove such results in a rigorous mathematical way we need to apply advanced methods 
from the geometry of nonholonomic manifolds. For our purposes, such manifolds can be considered as usual 
pseudo-Riemannian spacetimes, endowed with additional non-integrable distributions and frame structures. 

Both classes of metrics m and © can be characterized, respectively, by scaling factors a(t) and a(t). 
Let us suppose that we have found a cosmological solution of type © in a given theory of modified gravity 
and analyze how this metric can be formally diagonalized for deformations of a small real parameter e (when 
0 < £ <S 1). We can consider “homogeneous” approximations of type /13(f) ~ 1 + 6X3 (t),Wi(t) ~ £Wi(t) and 
rii ~ enjJl In explicit form, such metric, with small off-diagonal deformations on e and rescaling a(t) —* a(f), 
can be written as 


ds 2 = a 2 (t)[(e 1 ) 2 + (e 2 ) 2 ] + a 2 (f)[l + £X 3 (t)](dy 3 + ehidx 1 ) 2 + (dt + £Wi(t)dx 1 ) 2 ■ (7) 

See below how it is possible to construct subclasses of off-diagonal configurations in a f (R, ...) gravity where 
Y (1321) goes into A (l!© . and <l 2 = A _1 [<f> 2 | Y| + f d( <h 2 cA| Y|] (Ifill) results in f —>• f = R, an effective 
R°g = A8°jg which admits LC-solutions with zero torsion. We will be able to reproduce the ACDM model 
provided the metric © defines certain classes of solutions constructed for a corresponding effective action 
in GR, namely 

s = ±f 5 4 u^/\ eg Q/J |( e R - 2A + m C( e g a/3 , )). (8) 

In this action, the Ricci scalar e R = R(a,e) is constructed for £ g a /3 with coefficients of ©, A is an effective 
cosmological constant used for nonholonomic deformations, and m C is considered for certain effective matter 
fields with certain pressure m p and energy-density m p. The EoS are chosen, for simplicity, to correspond 
to an effective de Sitter configuration determined by A, where w := Pa/pa = — 1, with pressure pA and 
energy-density pA- 

We can describe the theories determined by © and © with respect to nonholonomic (non-integrable) 
dual frames e“ = (e*,e a ), which is convenient for constructing off-diagonal solutions, or to redefine the 
constructions with respect to local coordinate coframes du a = (dx l ,dy a ), where certain analog of the FLRW 
metric and ACDM like theories can be analyzed. For £ —>• 0, the metric © transforms into 

ds 2 = a 2 (/)[(e 1 ) 2 + (e 2 ) 2 + (dy 3 )} 2 + dt 2 . (9) 

5 On inhomogeneity effects in cosmology, see |48l . In a more general context, it is possible to consider also “small” local 
anisotropic deformations depending on space like coordinates when \ 3 {x k , t), Wi ( t ) ~ ewi(x k , t) and rii ~ ehi (x k ). Some amount 
of anisotropy is compatible with observational data in various gravity and cosmological theories. See |491 !50j , for reviews of 
various approaches related to GR and generalizations of Bianchi, Kasner and Godel type configurations; [44l [45 ; 46 ; 47] . for 
off-diagonal configurations; and [51 j . for f{R )-modified gravity theories. We note also that the approximation / 13 (f) « 1 + EX 3 R) 
can be very restrictive—one can consider more general classes of solutions with arbitrary /i 3 (t). 
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which is just ID but with a re-scaled factor because of the nonholonomic transformations T —>• A and 

$ ->• < i >. 

The corresponding Einstein equations with respect to the nonholonomic frames are 

2>H 2 = k 2 m p + A, (10) 

2 H° = -k“( m p + m P + A), 

where H° := a°/a. We can express £ R + m £= a R+ ^C+e ^C , where a R and are computed for the 
metric © and include all e-deformations in (f8j) . Such } n C results in the effective splitting m p = ^ p+e ^p 
and m p = ® n p + e ^p. In this way, we can encode the off-diagonal components as certain additional terms 
into the matter source, or either consider them as a polarization of the effective cosmological constant 
A := A + £ 1 A@ We also note that possible small inhomogeneous and locally anisotropic contributions, and 
concordance with observational data, can be estimated similarly to those presented, e.g., in [51] . This could 
be a ground for further investigations of such “slightly” /(-R)-modified off-diagonal cosmological models. 
In this subsection we provide only a few qualitative estimations, in order to demonstrate that “realistic” 
ACDM like cosmological theories can be equivalently modeled both in terms of nonholononomic frames and 
of coordinate frames, if small off-diagonal deformations are considered, only, and then the limit e —>• 0 is 
taken. 

In coordinate frames, Eqs. (HOD are written as 

3 H 2 = K 2 ° m p + A, 

2 H« = -K 2 (° n p+ ° m P + A). 

For e —> 0, the diagonalized solutions are determined by a (and not by a in ©) and can be parameterized to 
define and effective ACDM like model where a = a c e Hct , for a positive constant a c . Thus, modified gravities 
with equivalent off-diagonal encodings of /(i?)-modified gravity seem to result in realistic cosmological 
models, at least for small parametric e-deformations. 

The main goal of this work is to study possible nonlinear gravitational and matter field interactions 
which result in the encoding of modified gravities into generic off-diagonal metrics defining effective Ein¬ 
stein spaces, without certain special assumptions on the linearization of some associated systems of PDEs 
and their solutions. Surprisingly enough, the AFDM allows us to find such “non-perturbative” solutions 
in explicit form, by using geometrical methods. The values a(t), and ni are nonlinearly deter¬ 

mined by generating functions and sources of type <F(i) and T(t). In general, such nonlinear modifications 
of a “prime” a(t) are not small. Even if we can introduce an effective scaling factor a(t), this value de¬ 
scribes a nonlinear and inhomogeneous evolution with respect to nonholonomic (nonintegrable) dual frames 
e“ = (e\e a ). All generic off-diagonal cosmological models can be also redefined with respect to local 
coordinate coframes du a = (dx l ,dy a ). In local coordinate form, we are not able to analyze common and 
different properties of diagonalizable and non-diagonalizable models only by comparing the evolutions of 
a(t) and a(t). From the physical point of view, we can consider the Universe as an aether with a complex 
vacuum and nonvacuum nonlinear structure determined by possible /(-R)-modifications. An observer ac¬ 
quires experimental/observational data with respect to a local comoving frame e Q = (e*,e a ) where generic 
off-diagonal gravitational and matter field interactions are taken into consideration. For certain parametric 
resonant dependencies, even the smallest nonlinearities can result in substantial polarizations of the gravita¬ 
tional vacuum aether, with possible Lie group or solitonic symmetries, or without any anisotropic symmetry 
prescribed in advance. Such cosmological models are described by more sophisticate geometries than the 
FLRW cosmology (for references, see |481 (50]). 

The key idea of our work is that, within certain assumptions, various possible /(i?)-nonlinear modifica¬ 
tions can be encoded into off-diagonal terms and some effective a(t), hs(t),Wi(t ) via nonlinear interactions. 

f, We do not provide here explicit formulas for the corrections proportional to e because, in the end, we shall take smooth 
limits e —» 0. The main constructions for nonholonomic off-diagonal transforms are based on rescaling a(t) —» a 2 (t) generated 
by the solutions with T —> A and $—>■$. 
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This can be done for more general classes of cosmological solutions with nonlinear gravitational interactions 
restructuring the spacetime aether before considering certain small e-parameters. Such nonlinear cosmologi¬ 
cal evolution is determined by three functions of a time like variable, t, characterizing a more complex model 
then the FLRW one. We get, indeed: 1) a scaling factor a(t): 2) a diagonal inhomogeneity function / 13 (f); 
and 3) off-diagonal deformations via Wi(t ). Having constructed a class of off-diagonal solutions then, with 
certain additional assumptions for the effective linearization in terms of e, one can study possible observable 
e-small inhomogeneous or locally anisotropic contributions. The physical effects of small e-deformations can 
be compared, for instance, with those for a scaling factor a(t), although this will not be the aim of this paper. 
The main results and conclusions of it (see Sects. 4 and 5) will have to do with certain special properties 
of off-diagonal nonlinear systems with nonintegrable constraints and with the exact solutions. Even for the 
physical interpretations of observable cosmological data at a fixed time t = to, we can take the limit e —>• 0 , 
where / 13 (f) —» 1 , and Wi and m may vanish or result in a nonholonomic frame structure; a generalized 
nonlinear cosmological evolution by such generalized solutions may result in a modified scaling factor a(t) 
encoding both /(./^-modifications and off-diagonal nonlinear interactions for t < to- 

2.2 Conventions and geometric preliminaries 
2.2.1 Nonlinear connections and N-adapted frames 

Let us consider a pseudo-Riemannian manifold V, dim V = n + m, (n,m > 2). A Whitney sum N is 
defined for its tangent space TV, 

N: TV = hTV © vTV. (11) 

Conventionally, this states a nonholonomic (equivalently, non-integrable, or anholonomic) horizontal (h) and 
vertical (v) splitting, or a nonlinear connection (N-connection) structure. I 11 local form, it is determined by 
its coefficients N = {N?(u)}, when N = N?(x,y)dx l < 8 > d/dy a for certain local coordinates u = (x,y), or 
u a = ( x i , y a ), and h-indices i,j, ... = 1, 2, ...n and u-indices a, b, ... = n + 1, n + 2,..., n + mQ Such a h-v- 
decomposition can be naturally associated with some N-adapted frame or, respectively, dual frame structures, 
e u = (e;,e a ) and = (e\e a ), 


e* = d/dx 1 - N?{u)d/dy a , e a = 8 a = d/dy a 
and e* = dx l , e a = dy a + N^(u)dx‘ l . 


( 12 ) 

(13) 


The nonholonomy relations hold 


[ e a, 6/s] ^^a/3^' Y> 


(14) 


with nontrivial anholonomy coefficients Wf a = d a N^, W°' t = = ej (ZV“) — ej(/V a ). The coefficients 


define the N-connection curvature. 


2.2.2 Distinguished metric structures 


Any metric structure g on V (for physical applications, we consider pseudo-Euclidean signatures of type 
(+, +, +, —)) can be written in two equivalent ways: 1) with respect to a dual local coordinate basis, 


g = g „du a <g) du ^, 

0 z-ap 


(15) 


where 


Qa/3 


9ij + N^N^gab N? g ae 
A) gbe 9ab 


(16) 


7 The Einstein rule on index summation will be applied if the contrary is not stated. For convenience, “primed” and “un¬ 
derlined” indices will be used, and boldface letters to emphasize that an N-connection spitting is considered on a spacetime 
manifold V = (V, N). 
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or 2) as a distinguished metric (in brief, d-metric, i.e. in N-adapted form, 


g = 9a(u)e a ® = gi(x k )dx l ® dx l + g a (x k ,y b )e a ® e a . (17) 

To prove the decoupling of fundamental gravitational equations in modified gravity is possible for d-metrics 
and working with respect to N-adapted frames. 


2.2.3 Distinguished connections 

A linear connection is called distinguished, d-connection, D = ( hD,vD ), if it preserves under parallelism 
a prescribed N-connection splitting ED- Any D defines an operator of covariant derivation, DxY, for a 
d-vector field Y in the direction of a d-vector X. We note that any vector Y(u ) € TV can be parameterized 
as a d-vector, Y = Y“e Q = Y*ej + Y a e a , or Y = (. hY,vY ), with hY = {Y*} and vY = {Y a }, where 
the N-adapted base vectors and duals, or covectors, are chosen in N-adapted form El) and (1 1 3|) . TIig 
local coefficients of DxY can be computed for D = {T 7 ^ = ,Cj c ,Cfc)} and h-v-components of 

D ea eg := D a e^ using X = e Q and Y = eglf| We can characterize a d-connection by three fundamental 
geometric objects: the d-torsion, T, the non-metricity, Q, and the d-curvature, 7 Z, respectively defined by 


T(X,Y) := D X Y — D y X — [X, Y], <2(X) := D x g, 
^■(X, Y) := D x Dy — DyD x — Dj x y]- 


The N-adapted coefficients, 


T _ r T-l7 _ (rpi rpi rpa rpa rpa \ \ Y) _ / /~v7 

’ — i- 1 - a/3 — \ J jki 1 jai 1 jii 1 bi’ 1 be) a ph 

_ fT3 a _ ( T?a ni pc p i pc \ i 

/V — ^5— hjki n bjk’ n hja' 1 ' bja’ n hbco n bea) h 


(18) 


of such fundamental geometric objects are computed by introducing X = e a and Y = e^, and D = {T 7 a ^} 
into the formulas above (see [271123 ES EH EH for details). 


2.2.4 Preferred d-metric and d-connection structures 

A d-connection D is compatible with a d-metric g if and only if Q = Dg = 0. Any metric structure 
g on V is characterized by a unique metric compatible and torsionless linear connection called the Levi- 
Civita (LC) connection, V. It should be noted that V is not a d-connection because it does not preserve 
under parallelism the N-connection splitting ED- Nevertheless, such a h-v decomposition allows us to define 
N-adapted distortions of any d-connection D, 


D = V + Z, 


(19) 


with respective conventional “non-boldface” and “boldface” symbols for the coefficients: V = {T 0 ^,} and, for 
the distortion d-tensor, Z = {ZV }. 

This stands for any prescribed N and g = hg + vg, but alternatively to V, on V, we can work with the 
so-called canonical d-connection, D, when 


(g<N) 


V : 

D : 


Dg 


Vg = 0; V T = 0; 

0; hT = 0, vT = 0, hvT / 0; 


are completely defined by the same metric structure. The canonical distortion d-tensor Z in the distortion 
relation of type fU, D = V + Z, is an algebraic combination of the coefficients of the corresponding 
torsion d-tensor T = {T"^}. The respective coefficients of the torsions, T and v T = 0, and curvatures, 

8 We shall use the terms d-vector, d-tensor, etc. for any vector, tensor valued with coefficients defined in a N-adapted form 
with respect to the necessary types of tensor products of bases, (El) and m , and necessary h-u-decompositions. 
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1Z = {R 0 ^} and v 7£ = {_R°jg of D and V can be defined and computed using formulas similar to (1181) . 
We note that the coefficients T°g are not trivial but nonholonomically induced by anholonomy coefficients 
W2p (1141) and certain off-diagonal coefficients of the metric (1161) FI 
The Ricci tensors of D and V are computed in the standard form, 

IZic = {R 

P'y R 7 a ,g7} an< ^ ^ C ~ P'r }' 


With respect to N-adapted coframes (1131) . the Ricci d-tensor IZic is characterized by four h-v N-adapted 
coefficients 

R«/3 = {Rij := R ijk -1 Ria := ~R ikai Rai := R aibi Rab := R afcc}) (20) 

and (an alternative to the LC-scalar curvature, R := g a ^R a p) scalar curvature, 


R:=g^R a p = g ij Ri j + g ab R ab - 


( 21 ) 


We emphasize that any (pseudo) Riemannian geometry can be equivalently described by both geometric 
data (g, V) and (g, N,D). For instance, there are canonical distortion relations 

1Z = v 1Z+ v Z and 7Zic = Ric + Zic, 

where the respective distortion d-tensors v Z and Zic are computed by introducing D = V + Z into the 
corresponding formulas (1181) and (1201) . The canonical data (g,N,D) provide an example of nonholonomic 
(pseudo-) Riemannian manifold which is a standard one but enabled with a nonholonomic distribution 
determined by (g, N). If the coefficients = 0, such a distribution is holonomic, i.e. integrable. 

Nevertheless, physical theories formulated in terms of data as (g, V), or (g,N,D), are not equivalent if 
certain additional conditions are not imposed. Let us consider an explicit example. We can introduce the 
Einstein d-tensor of D, 

1 


Ra/3 •— R»/3 2 §o/3 R, 

and construct a N-adapted energy momentum tensor for a Lagrange density 

™ 2 

J-a/3 •“- 


( 22 ) 

1 C, of the matter fields, 


-C) 


^ ; performing a N-adapted variational calculus with respect to N-elongated (co) 

frames m and (USD, and consider that D is used as covariant derivative instead of V. In this way a 
nonholonomic deformation of Einstein’s gravity is constructed, being V —>• D = V + Z, with gravitational 
field equations 

Rq/ 3 = r2 (T a/3 — 2 § c*/?T) (23) 

for a conventional gravitational constant k 2 and T := Such equations are different from the 

standard Einstein equations in GR because, in general, R Q( 3 ^ R a p and T a p ^ T a p, where T a a : = 

LC-configurations can be extracted from certain classes of solutions of Eqs. (123[) if additional conditions 
are imposed, resulting in zero values for the canonical d-torsion, T = 0. In N-adapted coefficient form, such 
condition is equivalent to 


rpi _ ji T i rpl _ r-pa _ (~\a rpC _ T C ( ]\TC\ rpa _ s~ia s^ia 

jk ^ jk ja ^ jb>ji ^ jii aj ^aj j ) > 1 be ^bc cb m 


(24) 


It should be emphasized that we are able to find generic off-diagonal solutions of the Einstein equations 
in GR depending on three and more coordinates for D —>• V, when R a /j —>• R a p and T a p —>• T a; g, if the 
nonholonomic constraints (1241) are imposed after certain classes of solutions were found for D^V. But we 
are not able to decouple such systems of nonlinear PDEs if the zero torsion condition for V is imposed from 
the very beginning. 

9 In the Riemann-Cartan geometry, such a torsion is for a general metric compatible linear connection, D , which is not 
necessarily completely defined by the data (g, N). 
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2.3 Nonholonomic structures in modified gravity theories 

In general, different models of modified gravity are formulated for independent metric and linear con¬ 
nection fields with a corresponding Palatini type variational formulation (see MM El mam 

mmm)- The gravitational and matter field equations in modified gravities consist in very sophisticate 
systems of nonlinear PDEs for which finding exact solutions is a very difficult technical task, even for the 
simplest, diagonal ansatze with the coefficients of the metrics and connections depending on just one (time 
or space variable). Nevertheless, the anholonomic frame deformation method [ 2 ?, [251 [29l l30l 1ST] seems to 
work efficiently and allows to construct off-diagonal solutions in modified gravity theories mmmm- 

2.3.1 Equivalent modeling of modified gravity 

Consider three classes of equivalent theories of modified gravity defined for the same metric field g = {g^} 
but with different actions (and related functionals) for gravity, g S , and matter, m S, fields, 

5 = *5+ m S = ^J f (R, T, R a pT a P) y/\g\d 4 u + J m CVW\d A u 
= g S+ m S = ^J ?(R,T,R Q/3 T^)^ljd 4 ir + J m L^\g\d 4 u 
= 9 S+ m S = ^f + kf (25) 

Here, we use boldface d 4 u in order to emphasize that the integration volume is for N-elongated partial 
derivatives (1131) . k 2 is the gravitational coupling constant, the values with are computed for a canonical 
d-connection D and the values with “V” for re-defined geometric data (g,N,D) for certain nonholonomic 
frame transforms and nonholonomic deformations g a $ ~ g a p ~ ga/30 For simplicity, we consider matter 
actions which only depend on the coefficients of a metric field and not on their derivatives, 

T aP = m L g a d + 2S( m L)/ Sgafi- 

Such variations can be performed with respect to coordinate frames for g^ u or in various N-adapted forms 
for g^ and g^. 

2.3.2 Off-diagonal deformations of FLRW metrics 

We assume that the matter content of the universe can be approximated by a perfect (pressureless) fluid, 
where 

T „/3 = PSa/3 + (P + p)v a vp (26) 

is defined for certain (effective) energy and pressure densities, respectively, v Q being the four-velocity of the 
fluid for which v Q v“ = —1 and v Q = (0,0,0,1) in N-adapted comoving frames/coordinates. Here frame 
deformations/transforms of metrics of type g a p = e a a e P pga'p'i w iH be studied, being the FLRW diagonalized 
element 

ds 2 = g a '/ 3 'du a du^ _ a 2 (t)[dr 2 + r 2 d6 2 + r 2 sin 2 Odtp 2 ] — dt 2 , (27) 

= d 2 (t)[dx 2 + dy 2 + dz 2 ] — dt 2 

where the scale factor a(t) (we use also the value H := a/a, for a? := da/dt ) with signature (+,+,+,—), 
and a parametrization of coordinates in the form u a = ( x 1 = r, x 2 = 6,y 3 = ip,y 4 = t), or as Cartesian 
coordinates {x l = x,x 2 = y, y 3 = z,y 4 = t). For such cosmological metrics, the main issues of the Einstein 
and modified Universes are encoded into energy-momentum tensor T a p = pg a p + (p+p)v a Vg (we omit primes 

10 We shall give details in Sects. I3ll.3l and l3.2l 
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or other distinctions in the coordinate indices if there is no ambiguity) arising from a matter Lagrangian m C 
through 

T{t) = T “ = -p, P(t) = R a pT a P = R AA T 44 = —3p(H 2 + H°), (28) 

for T a p = diag[0,0,0,—p\. 

We will consider nonhomogeneous and locally anisotropic cosmological solutions of type (11611 and/or (El) 
generated by off-diagonal deformations of ED 

9i = 9 i(x k ) = rii(x k ,y 4 )gi(x k ,y 4 ) = e^ xk \ (29) 

g a = u 2 (x k ,y 4 )h a (x k ,y 4 ) = u 2 (x k ,y 4 )i) a (x k ,y 4 )g a (x k ,y 4 ), 

Nf = rii(x k ),N 4 = Wi(x k ,y 4 ). 

In Eqs. (12911 there is no summation on repeated indices, g a = (pi,g a ) are polarization functions, the N- 
connection coefficients are determined by m and Wi, the vertical conformal factor uo may depend on all 
spacetime coordinates and g a = (gi,g a ) define the “prime” diagonal metric if rj a = 1 and iV“ = 0. The 
“target” off-diagonal metrics are with Killing symmetry on d/dy 3 when the coefficients (l29l) do not depend 
on 2/ 3 0 

ds 2 = a 2 (x k ,t)[rji(x k ,t)(dx 1 ) 2 + 772 (x k ,t)(dx 2 ) 2 ] + a 2 (x k , t)h^(x k , t)(e 3 ) 2 + uj 2 (x k , t)h A (x k , f)(e 4 ) 2 , (30) 

when a 2 (x k , t)rn(x k , t) = e^ xk \ for i = 1,2; a 2 h-j = uj 2 (x k , t)hs(x k , t), and 

e 3 = dy 3 + rii(x k )dx l ,e 4 = dy 4 + Wi(x k ,t)dx l . 

Functions gi,rj a ,a,ii>,uj,ni,Wi will be found such that, via nonholonomic transforms (12911 . when g a >p>(t) (I27|) 
—>• % a p{x k ,t) (1301) . off-diagonal nonhomogeneous cosmological solutions are generated in a model of modified 
gravity ED. We can consider subclasses of off-diagonal cosmological solutions but with deformed symmetries 
when certain nontrivial limits g a p(x k ,t) —>• g a pif) can be found and define viable cosmological models. 


2.3.3 Field equations for nonholonomic modified gravities and FLRW cosmology 

Applying an N-adapted variational procedure with respect to a nonholonomic basis (1 12[) and (1T31) for the 
action S = 9 S + m S, which is similar to that in [23, 23 EE] but for V —>■ D and matter source T a p (l26l) . 
we obtain the field equations for the corresponding modified gravity theory 


1 f — ^ ga/3f + (ga/SD^D/i — DqD^) 1 f + (T a 0 + ©a/?) 2 f + 

H a /3 3 f + I(D^T^ 3 f + g a /jD M D„T'“' 3 f) - D„D (a T^ 3 f = k 2 T 


(31) 


for ® a /3 p Sa/3 2T q/3 , 2 E p'jy p E a g 

with respective d-tensors defined by Eqs. EU1) . (Ell) and ED ; where 3 f := <9f/<9R, 2 f := di/dT and 
3 f := <9f/3P, when P = R Q/ 3 T Q ^ and (a/3) denotes symmetrization of the indices. 

In general, the divergence with D and/or V of Eqs. (1311) is not zero. Also Eqs. (l23l) have a similar 
property. In the last case, we can obtain the continuity equations as in GR and then deform them by using 
the distortions (1191) . which for the canonical d-connections are completely determined by the metric structure. 
There are certain types of conservation laws for matter fields with additional nonholonomic constraints. 
We can consider the field equations m and the equations derived by taking the divergence with D as 

11 We can consider nonholonomic deformations with non-Killing symmetries when, for instance, uj(x k ,y 4 ) —> uo(x k ,y 3 ,y 4 ), 
which results in a more cumbersome calculus and geometric techniques. For simplicity, we do not study such generalizations in 
this work (see examples in |2711281 29] 30] 130j 1. 
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nonholonomic distortions of similar systems of nonlinear functional PDEs considered in [24, [251 E5] for V. 
Remarkably, such sophisticate nonholonomic and nonlinear systems can be solved in very general off-diagonal 
forms, by applying the anholonomic frame deformation method. In order to compare these results and to 
find possible applications in modern cosmology, we will consider a particular equation of state (EoS) p = wp 
with w = const, and study the cosmology of off-diagonal distortions of certain FLRW models considered in 
the framework of GR and its modifications. In both cases, by exploring some particular classes of solutions, 
the dynamics of the matter sector of generalized f{R, T, R^ U T^ V ) gravity (with respect to N-adapted frames) 
may lead to similar cosmological scenarios as GR, but with nonholonomic constraints and deformations. 


3 The anholonomic frame deformation method and exact solutions in 
modfied gravities 

A surprising property of Eqs. (I23p and (1311) is that they can be integrated in very general form with 
generic off-diagonal metrics when their coefficients depend on all spacetime coordinates via various classes 
of generating and integration functions and constants. In particular, we can consider such generating and 
integration functions when g a p(x k ,t) (1301) result in off-diagonal metrics of type g a fi{t) depending on the 
parameters and possible (non-) commutative Lie algebra or algebroid symmetries. 


3.1 Off-diagonal FLRW like cosmological models 


We shall study cosmological models with sources of type fl26l) when the four-velocity v a is re-parameterized 
in a way that for some frame transforms as 


y a B '■= K 2 (T a p — -go^T) 

diag[T 1 = T 2 , T 2 = h r(x% T 3 = T 4 , T 4 = v T(x\t)\ (32) 

—>• A g a p (redefining the generating functions and sources), (33) 


for effective h- and u-polarized sources, respectively, h T(x i ) and T 4 = V T(x 2 ,t), or an effective cosmological 
constant A. For simplicity, we can consider effective matter sources and “prime” metrics with Killin g sy mmetry 
on d/ds, i.e. when the effective matter sources and d-metrics do not depend on the coordinate y 3 o In brief, 
the partial derivatives d a = d/du a on a 4-d manifold will be written as s* = ds/dx 1 ,s / = ds/dx 2 ,s* = 
ds/dy 3 , s° = ds/dy A . 

The nontrivial components of the Ricci d-tensor (1201) and nonholonomic Einstein equations (1221) . with 
source (1321) parameterized with respect to N-adapted bases (ED and (ED, for a d-metric ansatze ED with 
coefficients (1291) . are 




Rzk 




~R 2 2 = 


-Ri = 


n 3 „oo 

2 /i 4 k 


1 


2ffi52 

1 

2/I3/I4 


[»;•- 


9i 92 idl) 


[ht~ 


2ffi 2 52 

(htf 


2/i, 


+ 4 2 3) 2 h 4 


w k , u oo jht) 
2/13 [ 3 2/13 


MM 1 

2/14 


9i 


MMi 

2/14 J 

= 0, 


9 W 2 

252 


= "T, 


h% .dkh'i 
4/i 3 /13 


dkhj 

/14 


Ml! i = 

251 J 


MM 

2/13 


(34) 

(35) 

(36) 

(37) 


where are computed for ui = 1 and then the formulas are generalized for uj 7 ^ 1 via v-conformal transforms 
(see Refs. |241 [25l[26] for details), 

eiU = diU — Hi co* — WiU 0 = 0. (38) 

12 The method can be extended to account for y 3 dependence and non-Killing configurations (see [27, 28j (29] 130)1. In this 
paper the local coordinates and ansatze for d-metrics are parameterized in different forms than in previous works, what is more 
convenient for the study of cosmological models. 
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The d-torsion (1241) vanishes if the (Levi-Civita, LC) conditions LT = e a (AT),Cj b = 0, fl'K = 0, are satisfied 
for 

Wi = ( 8 i - Wid A ) In ^S\^ul (di - Wid 4 ) In \/\h^\ = 0, (39) 

d kWi = diW k ,rii = 0,5 in k = d k rn. 

The above system of equations can be integrated in very general situations, for instance, for d-metrics with 
Killing symmetry on 83. 


3.1.1 Decoupling of PDEs for inhomogeneous cosmological metrics 

The system of equations (134D - (1381) has an important decoupling property. To show this explicitly, we 
rewrite it as nonlinear PDE which posses an important decoupling property, allowing integration step by 
step of such equations. For 7 ^ 0, h T, V T 7 ^ 0, Killing symmetry on 83 and parameterizations (l29l) . these 
equations can be written as 


= 2 h T (40) 

(As = 2 /i 3 /i 4 "T (41) 

nV + in* = 0, (42) 

(3wi - at = 0, (43) 

diLd - (<9;0/</P)uP = 0, (44) 

for o 

oti = h%8i(t),j3 = h% 0 °,7 = (ln|/i 3 | 3/2 /|/i 4 |) , (45) 

where 

(j) = In |/i 3 /\/1/i 3 / 14 11, and/ or <h := e(46) 


is considered as a generating function. Eq. (1441) is just Eq. (l38l) for a nontrivial solution of (l43l) with 
coefficients (145[) . when 


Wi = di^/cj) 0 . 


(47) 


The decoupling property of the above system of equations follows from the facts that: 1) integrating the 2-d 
Laplace equation (140[) one finds solutions for the /i-coefficients of the d-metric, and 2) the solutions for the 
coefficients of the d-metric can be found from (1411) and (1461) . 3) Then the N-connection coefficients Wi and 
Hi can be found from (1421) and (1431) . respectively. 


3.1.2 Cosmological solutions with nonholonomically induced torsion 

The equations (1401) and (143[) can be solved, respectively, for any source h T(x k ) and generating function 
(j)(x k 1 t). The system (1411) and (l46l) can be written under the form 

h 3 h 4 = <As /2 "T and \h 3 h A \ = (^)V 2 ^, 


for any nontrivial source 
(e 2 ^)® _ „ u 

4| 2 1 v T| > 1 ' e ‘ '*'3 

equation, we obtain 


'Tin (1411) . Introducing the first equation into the second, one finds |/i|| = 
°/r 3 (x fc ) + f . where °/i 3 (x fc ) and € 3,64 = ±1. Using again the first 


h>A = 


</°(ln y / M) C 
2 


1 




2 B T $ h -3 


(48) 


We can simplify such formulas for /i 3 and /14 if we redefine the generating function, —>• <F, where 

($ 2 )°/l « T | = ($2)o/ A> j e 


<F = A 


-1 


$ 2 


T| + dt T 2 | u T| 


(49) 
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for an effective cosmological constant A which may take positive or negative values. We can integrate on t, 
include the integration function °/i3(x ,fc ) in $ and write 

/i 3 [S] = $ 2 /4A. (50) 


Introducing this formula and (1491) into (1481) . we compute 

(ln|$|)° (8 2 )° 


M$] = 


4| «T| 


$ 2 | W T| + J 


i -i 


dt $ 2 | 


As next step, we need solve Eq. (I42[) by integrating on t twice. We obtain 

n k = l n k + 2 n k J dt h A /(^/\h^\) 3 , 


(51) 


(52) 


where i n k (x i ), 2 'n k {x' 1 ) are integration functions and h a [tf>] are given by formulas (1501) and (1511) . If we fix 
2 n k = 0, we shall be able to find n k = i n k {x l ) which have zero torsion limits (see examples in subsection 

MM - 

The solutions of (l43l) are given by (1471) , which for different types of generating functions are parameterized 


as 


dj$ _ dj{§ 2 ) 

$0 “ ($ 2 )* ’ 


(53) 


where the integral functional < h[<I > , v T] is given by (l49l) . 

We can introduce certain polarization functions rj a in order to write the d-metric of such solutions in 
the form (1301) . Let us fix u 2 = | /j-41 —1 to satisfy the condition (1441) . which for a generating function $[0] is 
equivalent to 

dih 4 - di<5> h\ = 0. (54) 


These first order PDE equations impose certain conditions on the class of generating function $ and source 
V T. For instance, we can choose such a system of coordinates where = -j(e~^)° which transforms Eq. ED 
into /14 = <h, i.e. this coefficient of the d-metric is considered as a generating function, and Eqs. (l54l) are 
solved. In general, the integral varieties of such equations cannot be expressed in explicit holonomic form. 

A modification of the scale factor a(t) a(x k ,t), for the FLRW metric (l27l) (with for g\ = g 2 = 
(?3 = a 2 ,g 4 = —1, has to be chosen in order to explain observational cosmological data. For any prescribed 
functions a(x k ,t ) and ui 2 = |1 — 1 and solutions \ (see (1401) 1 and h a [4>], n k (x l ), u;j[d ) ] (given respectively 
by formulas (l50l) - (l53j) l. we can compute the polarization functions rji = a~ 2 e^,r ]3 = d~ 2 hz,r]i = 1 and 
function /13 = ^.3/a. 2 1 /x4[. Such coefficients (see the data (1291) 1. define off-diagonal metrics of type (I30|) . 


ds 2 = ^(x^ ,t)[gi(x lc ,t)(dx 1 ) 2 + rj 2 (x lc ,t)(dx z ) z } + a z (x >: , t)h 3(x fc , t)[dy 6 + rii(x k )dx l ] z — [dt + 


„ 2\2 


il2 


dM$, V T] 

«T] 


dx 


i i2 


. _ (55) 

Choosing any generating functions a 2 (x k , t), ip(x i ) and < L[ ( I > , U Y] and integration functions ni{x k ), we generate 
a nonhomogeneous cosmological model with nonholonomically induced torsion (|241) . More general torsions 
can be induced if rii(x k ,t ) is taken with two types of integration functions 1 rii(x k ) and 2 ni(x k ) (see 
Eqs. (1521) 1. Having constructed this solution, we can now consider certain subclasses of generating and 
integration functions where a(x k ,t ) —» a(t) 7^ a(t),Wi —>• Wi(t),rii —>• const, etc. In this way generic off- 
diagonal cosmological metrics are generated (because there are nontrivial anholonomy coefficients W^ a in 


©). 
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3.1.3 Levi-Civita off-diagonal cosmological configurations 

The LC-conditions (1391) are given by a set of nonholonomic constraints which cannot be solved in explicit 
form for arbitrary data (<f>,T) and integration functions iand 2 r ifc- However, some subclasses of off- 
diagonal solutions can still be constructed where via frame and coordinate transforms we can chose = 0 
and infc = d^n with a function n = n(x k ). It should be noted that (di — Wid^Q = 0 for any &(x k , y A ) if Wi 
is defined by (1531) . Introducing a new functional -£>(<!>), we find that (di — Wjd^)B = (dj — Wid^jQ = 0. 
Using Eq. (l50l) for functionals of type h 3 = .B(|5>(3>)|), we solve Eqs. (dj — Widths = 0, what is equivalent 
to the second system of equations in (l39l) . because (di — w^d 4 ) In y/\hz f (di - Wid 4 )h 3 . 

We can use a subclass of generating functions <f> = $ for which 




(56) 


and get for the left part of the second equation in (1391) . (di — W{d 4 ) In \/\h 3 \ = 0. The first system of equations 
in (1391) can be solved in explicit form if Wi are determined by formulas (1531) . and h 3 [<&] and /i 4 [<I>, <L 0 ] are 
chosen respectively for any T —>• A. We can consider <f> = 4>(ln ^/|h- 4 |) for a functional dependence /i 4 [$[$]]. 
This allows us to obtain Wi = <9j| < h|/|<l >| 0 = di\ In 1/i 4 11/1 In y^/^H 0 . Taking the derivative 84 on both sides 
of these equations, we get 

O _ (di\ In VWillT I In V«r 

W *~ IlnVTMl^ "|lnv«r- 

If the mentioned conditions are satisfied, we can construct in explicit form generic off-diagonal configurations 
with w i = (di — Wid 4 ) In ^/jhlj", which is necessary for the zero torsion conditions. 

We need als solve for the conditions d^Wi = diW from the second line in (I39[) . We find in explicit form 
the solutions for such coefficients if 

Wi = di = diA, (57) 

with a nontrivial function A(x k ,y 4 ) depending functionally on the generating function T. 

Finally, we conclude that we generate LC-configurations for a class of off-diagonal cosmological metric 
type (1171) for T = T = A, 4> = <L = and 2 = 0 in (l52l) which are parameterized by quadratic elements 


ds 2 = e^[(dx 1 ) 2 + (dx 2 ) 2 ] + 


^ [ dy 3 + (5 fc n(* i ))ds fc ] 2 - ^ 


4|A| 


| A|ff > 2 


[dt + (diA[Q])dx 


i i 2 


(58) 


We can re-write such solutions in the form (155)) . which provides us a general procedure of off-diagonal 
deformations with a(t ) —>• a(x k ,t ) (see the FLRW metric (|27jl ). resulting in nonhomogeneous cosmological 
metrics in GR. Prescribing a function a(x k ,t), a generating function <h(a: fc , t) satisfying the condition (1561) and 
a solution e^ x ) (see (1401) ), we respectively compute the v-conformal factor and the polarization functions 
for 


h 3 = h 3 /a 2 \h 4 \ = T 4 /4a 2 (<f>°) 2 , w 2 = |h 4 | 1 = |A|$ 2 /(I> 0 ) 2 , ry = a 2 e v ’,r ? 3 = a 2 h 3 = <f> 2 /4|A|a 2 ,774 = 1. 
Such coefficients (see data (l29l) ) transform the off-diagonal cosmological solutions (l58l) into metrics of type 

m, 

ds 2 = a 2 (x k ,t){[r] 1 (x k ,t)(dx 1 ) 2 + 7 ] 2 (x k ,t)(dx 2 ) 2 ] + h 3 (x k ,t,)[dy 3 + (dkn(x 1 ))dx k ] 2 }-[dt+(d i A[®])dx 1 } 2 . (59) 


The dependence on the source A is contained in explicit form in the polarization 773 , for instance. This class 
of effective Einstein off-diagonal metrics g a / 3 (x k ,t) define new nonhomogeneous cosmological solutions in 
GR as off-diagonal deformations of the FLRW cosmology. For certain well-defined conditions, one can find 
limits g a y —>• g a y(t,a(t),h 3 (t),^(t),yi(t)). This provides explicit geometric models of nonlinear off-diagonal 
anisotropic cosmological evolution which, with respect to N-adapted frames, describe a(t) with modified 
re-scaling factors. 


16 






3.2 Effective FLRW cosmology for /-modified gravity 

The anholonomic frame deformation method outlined in previous subsections can be applied for the 
generation of off-diagonal cosmological solutions of field equations of modified gravities, see (EH). Redefining 
the generating functions via the transforms (149 D and (I56p . —» <f> —>• <f>, we can generate off-diagonal 
cosmological configurations with R = 4A, see (1321) and (1331) . Such parameterizations of geometric data 
and sources are possible for certain general conditions via transforms of N-adapted frames when the action 
functional functionally depends on A and on the effective sources, f[R(A), T(A), P], with P(t) = R a ^T“^ = 
—3 p(H 2 + H°) and H = a°/a with scaling factor a(t ) taken for some limits of a solution (I Hop , or (I59p . 
It should be noted that in the variables corresponding to the Levi-Civita connection V the functional f —>• 
f(R,T,R^ u T^ v ) describe very general modifications of GR which in our approach are encoded into a very 
sophisticated off-diagonal effective vacuum structure with nontrivial vacuum constants. 

We assume that the density of matter p = p in T a p (l26l) is the same as for a standard FLRW metric 
m and does not change under off-diagonal deformations with respect to N-adapted frames. For such 
configurations, the functions and S Qj g are parameterized, respectively, as 

0/ = (P - 2A )5% 3/3 = (2A 2 —p A — ±4A 2 )^ = - P A6% 

where terms with A 2 compensate each other in 4-d. We can write D^T Q( g = 0, 1 f ~ <9 2 f/(9R 2 e^A ~ 0, 
and (similarly) 2 f ~ 0, 3 f ~ 0, for ~ T Q( g ~ A d a p, A = const, with respect to corresponding 

classes of N-adapted frames. Eqs. (1311) transform into a system of nonholonomic nonlinear PDEs of type 
El-Ell), r% = T<5 a p, with effective diagonalized source 

A i? 2 f 3 f 

T = -A + —- + (2A-n- 2 A-p)-k+pA— , (60) 

r f H >f 

which can be parameterized with dependencies on ( x l ,t ), or on t. These equations can be solved for very 
general off-diagonal forms, depending on generating and integration functions, following the procedure out¬ 
lined in previous subsections. Redefining the generation function as in (l49l) . when an effective cosmological 
constant A is generated from Tone has 


$ 2 = A 


-l 


$ 2 | XI + dt $ 2 | XT 


(61) 


Such a generating function defines off-diagonal cosmological solutions of type (1 55 P . or (l58l) . as solutions of 
held equations for an effective (nonholonomic) Einstein space RA = A8 a g. In this way, a geometric method is 

provided when the (effective or modified) matter sources transform as X (l32l) — >• A (l33l) and the gravitational 
held equations in modified gravity can be effectively expressed as nonholonomic Einstein spaces when the 
d-metric coefficients encode the contributions of f, 1 f, 2 f and 3 f and of the matter sources. 

We can consider inverse transforms with A —>• X for (1611) and state that for certain well-defined conditions 
of type (1561) and (1571) we can mimic both /-functional contributions and/or massive gravitational theories 
mi no mi nz]. Here we emphasize that off-diagonal configurations (of vacuum and non-vacuum types) are 
possible even if the effective sources from modihed gravity are constrained to be zero. 


4 Off-diagonal modeling of cosmological modified gravity theories 

This section has three goals. The hrst is to provide a reconstruction procedure for off-diagonal effective 
Einstein and modihed gravity cosmological scenarios. The second is to apply these methods in practice and 
provide explicit examples related to f(R) gravity and cosmology. The third goal is to analyze how matter 
stability problems for f(R )~theories can be solved by nonholonomic frame transforms and deformations and 
imposing non-integrable constraints. 
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4.1 Reconstructing nonholonomic general f(R )~models 


Let us construct an effective Einstein space which models a quite general modified gravity theory with 
f(R,T,R a pT a P) = R + F{R a pT a P) + G(T). This theory admits a reconstruction procedure which does 
not affect the observational constraints when a realistic evolution is studied |241 [25l [26 ]. Following the 
anholonomic frame deformation method with an auxiliary canonical d-connection D, the modified gravity 
m is formulated for 

f(R,T,R a/3 T a/3 ) = R + F(P) + G(T). (62) 


We can self-consistently embed this model into a nonholonomic background determined by N-adapted frames 
ca and (1131) for a generic off-diagonal solution (l59l) with limits D — » V and g a p —» g Q , j g(t, a(t), /13(f), <t(t),r]i(t)). 
With respect to such frames, the nonholonomic FLRW equations are similar to those found in section III B 

°f [23] 25, 26] (see the second paper for details on methods of constructing solutions and speculations on the 

[13 


problem of matter instability) 


4.1.1 /-modified off-diagonal FLRW equations 


The effective function a(t) defines in our case off-diagonal cosmological evolution scenarios which are 
different from those where a(t) stands for a standard diagonal FLRW cosmology. For H := a°/a, 1 G := 
dG/dT and 1 F := dF/dP , we have 


3 H 2 + l 


f + G - 3(3R 2 - FT) p X F - p(k 2 - X G) = 0, 


(63) 


^ \ 00 / ^ \ O ^ 

—3H 2 — 2H° -[f + G— (p 1 Fj -4H(p l F) - [3H 2 + FT) p *F] = 0 


An observer is here in a nonholonomic basis determined by 1V“ = {rii,Wi(t)} for a nontrivial off-diagonal 
vacuum with effective polarizations p a (t), and can test cosmological scenarios in terms of the redshift l + z = 
a~ 1 (t ) for P = P(z) and T = T(z), with a new “shift” derivative when (for instance, for a function s(t)) 
s* = -(1 + z)Hd z . 

The system of two equations (l63l) simplifies by extending it to a set of three equations for four unknown 
functions {f (z), G(z), p(z), <^(z)} with a new variable <j(z) := p X F, 

3H 2 + i[f (z) + G(z)] - |[3 H 2 - (1 + z)H(d z H)\ q(z)^H 2 ( 1 + z)d z q{z) - n 2 p(z) = 0, 

-3 H 2 + (1 + z)H(d z H) - l{f (z) + G(z) - [3 H 2 - (1 + z)H(d z H)\s(z) 

+ [3(1 + z)H 2 - (1 + z)H(d z H)\d z <;{z) + (1 + z) 2 d 2 zz <,(z)j = 0, 

(d z 1 F) s(z) - p(z) (d z f) = 0. (64) 

Here, by re-scaling the generating function, we have fixed the condition d z 1 G (z) = 0. Such a nontrivial 
term must be considered if one wants to transform f into a standard theory f(R, T, R a pT a P). The functional 
G(T), in both holonomic and nonholonomic forms, encodes a new degree of freedom for the evolution of the 
energy-density of type 

P = /+a- 3(1+ro) = Po(l + z)a 3 ^\ (65) 

which is taken for the dust matter approximation w when the evolution reduces to p r\_/ (1 + z) 3 . For the 
assumption that such an evolution can be considered with respect to N-adapted frames, the solutions of (I64[) 
are determined by data {f (z), G(z),?(z)} by replacing the second and third equations into the first one and 
obtaining a single fourth-order equation for f(z). 

13 In section III A of that work, a model with G(T) = 0 was investigated in detail. The conclusion was that in order to 
elaborate a realistic evolution it is necessary to consider nontrivial values for G(T). In nonholonomic variables, such term G(T) 
allows to encode f(R) modified theories and related into certain off-diagonal configurations in GR, which simplifies the solution 
of the problem of matter instability (see subsection 14.31 . 
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4.1.2 Reconstructing f-models and effective Einstein spaces 

The reconstruction procedure is restricted to fluids without pressure when such approximation is consid¬ 
ered locally with N-adapted frames and the expressions (|28D for (a, H , p) are re-defined in terms of (a, H , p)\ 
data are written with a script “0” if z = zq, with £ = k 2 po/3Hq. One should not confused, e.g., H and 
H 0 , because these values are computed for different FLRW solutions, with a{z) determined for a diagonal 
configuration and a(z) for an off-diagonal one, respectively. We can express 

O 7T2 orr2 

T = T“ = -^^(1 + zf and P = R Q g T Q/3 = — 3^—^(l + zf[H 2 - (1 + z)H(d z H)\. 

Following the approach outlined in Sect. IIIB of |24l 1251 [26] , we introduce the parameterizations 

F(P) = R 0 2 F(P) and G(T) = R 0 2 G(T), ( 66 ) 

where P =P/Pq and T = T/Tq, for Pq = —9iLg£/At 2 and To = —3 Hq£/k 2 . In correspondingly N-adapted 
variables, the off-diagonal cosmological solutions can be associated with a class of de Sitter (dS) solutions 
with effective cosmological constant A (see (1611) 1. where H(z) = Hq results in P = T = (1 + z ) 3 for the 
energy-density (1651) . In these variables, the solutions of (1641) can be written as 

F = ciP 61 +P {, 2 / 3 [c 2 cos(— lnP) + c 3 sin(— InP)] + c 4 + 3£P, (67) 

3 3 

G = ciT 61 +T 62 / 3 [c 2 cos(— hiT) +c 3 sin(^ In P)] + c 4 -3£T, 

3 3 

being the constants b\ = — 1.327,6 2 = 3.414 and 63 = 1.38. The values ci,c 2 ,C 3 and C 4 are integration con¬ 
stants, and the second set of constants ci,c 2 ,C 3 and C 4 can be expressed via such integration constants, and 
&i, & 2 and 63 . We omit explicit formulas because for general solutions they can be included in certain gener¬ 
ating or integration functions for the modified gravity equations and ultimately related to real observation 
data for the associated cosmological models. 

For off-diagonal configurations, the f(R, T, R a ^T“^) gravity positively allows for dS solutions in presence 
of non-constant fluids, not only due to the term P = Rq,^T“^ in (1251) . and respective gravitational field and 
cosmological equations. This is possible also because of the off-diagonal nonlinear gravitational interactions 
in the effective gravitational models. It should be emphasized that the reconstruction procedure elaborated 
in [231. [25l 126] . see also references therein, can be extended to more general classes of modified gravity theories, 
to Finsler like theories and the ensuing cosmological models Busang. Introducing (1671) and (1661) into (1621) . 
we reconstruct a function f = R + F(P) + G(T). As a result, we can associate an effective matter source T 
(I60p . which allows the definition of a corresponding generating function (16 ip (see also and (I49p l. Finally, 
we can reconstruct an off-diagonal cosmological solution with nonholonomically induced torsion of type (1551) 
or to model a similar cosmological metric for LC configurations (1551) (equivalently, (l59l) h 

4.2 How specific f(R) gravities and the FLRW cosmology are encoded in nonholonomic 
deformations? 

It is well known that any FLRW cosmology can be realized in a specific f(R) gravity (see Ref. [52l, 53j 
and, for further generalizations, |2411251 l26] d 14 l In this subsection we analyze two examples of reconstruction 
of /(R)-gravities where the "e-folding” variable £ := lna/oo = — ln(l + z) is used instead of the cosmological 
time t and in related nonholonomic off-diagonal deformations. For such models, we consider f = f(R) in 
(I25p . use Y(x*,£) = A/ + f/2 1 f instead of (I60p and introduce these values in Eq. (1611) . which can be 
parameterized with dependencies on (x*,C) (in particular, only on £), T 2 = A - 1 [<h 2 | Y| + f d( <h 2 c^| Y|], 
when d c * = d/dQ with s° = Hd^s for any function s. The matter energy density p is taken as in (1641) . 

14 We use a system of notations different from that article; here, e.g., N in used for the N-connection and we work with 
nonholonomic geometric objects. 
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We restrict ourselves to N-adapted frames, m and (USD, determined by an off-diagonal cosmological 
solution of the (modified) gravitational field equations, and can repeat all computations leading to Eqs. (2)- 
(7) in |52l 53] and prove that a modified gravity with f(R) realizes the FLRW cosmological model. Such 
solutions depend on the above source type Y(&*,£) and generating function d>(x l ,£); also the nonholonomic 
background can be modeled to be nonhomogeneous (via Wi and rq depending respectively on x 1 and £, or 
only on £). The off-diagonal analog of the field equation corresponding to the first FLRW equation is 

f(R) = (H 2 + H S C R)<9 C [f(R)] - 3 6H 2 [4 H + ( d c H ) 2 + Hd^H] d? c [f(R)]+K 2 p. 


In terms of an effective quadratic Hubble rate, q(C) '■= R 2 (C), and considering that £ = £(R) for certain 
parameterizations, this equation yields 


f(R) = -18 g (C(R))[d c 2 cg (C(R)) + 4%(C(R))] 


d 2 f(R) „ r 1„ df(R) 

—+ 6 fe(c( R )) + Man))}^- 

r/R- * cm. 


2 poa { 


— 3(l+tar)^ —3(l+ro)C(R) 
( 68 ) 


We can construct an off-diagonal cosmological model with metrics of type (1551) and nonholonomically induced 
torsion (when t —>• £) if a solution f(R) is used for computing Y and <!>. Modeling such nonlinear systems 
we can consider solutions of the field equations for an effective (nonholonomic) Einstein space R °/J = 

when certain terms of type df(R)/dR and higher derivatives vanish for a functional dependence f(A) with 
<9<-A = 0. The nonholonomic cosmological evolution is determined by off-diagonal coefficients of the metrics 
and by certain non-explicit relations for the functionals variables, like c/(£(R(A))) and (effective/modified) 
matter sources transform as Y (1321) —> A (1331) . 

LC-configurations can be modeled by off-diagonal cosmological metrics of type (1581) when the zero torsion 
conditions (1391) are satisfied. We obtain a standard expression (see [52l 153]) for the curvature of V, 


R = 3 %(£) + 12g(C), (69) 

if the polarization or generating functions for (1581) and the solutions of (1681) are taken for diagonal configu¬ 
rations. 


4.2.1 Off-diagonal encoding of f(R) gravity and reproduction of the ACDM-era 

We here provide an example of reconstruction of models of f(R) gravity and nonholonomically deformed 
GR when both reproduce the ACDM era. For simplicity, we do not consider a real matter source (if 
such a source exists, it can be easily encoded into a nontrivial vacuum structure with generic off-diagonal 
contributions). 

With respect to correspondingly N-adapted frames and for a(£) and H {£) determined by an off-diagonal 
solution (]55[l . with nonholonomically induced torsion, or (1591) . for LC-configurations, the FLRW equation for 
ACDM cosmology is given by 

3 k~ 2 H 2 = 3k~ 2 Hq + poa~ 3 = 3k~ 2 Hq + po% 3 e -3 ^. (70) 

This equation looks similar to the one for Einstein gravity for diagonal configurations but contains values 
determined, in general, for other classes of models with off-diagonal interactions. Thus, Hq and po are fixed to 
be certain constant values, after the coefficients of off-diagonal solutions are found, and for an approximation 
were the dependencies on (x l ,() are changed into dependencies on £ (via a corresponding re-definition of 
the generating functions and the effective sources). We can relate the first term on the rhs to an effective 
cosmological constant A (j33[) . which in our approach appears via a re-definition (1491) . The second term in 
the formula describes, in general, an inhomogeneous distribution of cold dark mater (CDM) with respect to 
N-adapted frames. In order to keep the similarity with the diagonalizable cosmological models in GR we can 
choose these integration constants for A = 12 Rq to survive in the limit Wi,rii —>• 0. It should be noted that 
such limit must be computed for “nonlinear” nonholonomic constraints via generating functions and effective 
sources. 


20 




Using (170D . the effective quadratic Hubble rate and the modified scalar curvature, R, are computed to 
be, respectively, 

q{ C) := H'o + At 2 p 0 «o 3 e^ 3< ’ and R = 3d^q(() + 12q(Q = 12 Hq + n 2 poa^' 3 e -3 ^. 

These functional formulas can be used for the dependencies on R if a necessary re-definition of the generation 
functions, or an approximation (x l ,() —>• C i s performed. Expressing 

3C = — ln[fiU 2 /9Q 1 Oq(R — 12Rq)] and X := —3 + R/3Rg, 

we obtain from Eq. (16811 

d 2 ? df - 

X(l - x )-^2 + [*3 - (xi + X2 + 1)-A] — - XiX2f = o, (71) 

for certain constants, for which xi +X 2 = X 1 X 2 = —1/6 and X ‘3 = —1/2- The solutions of this equation with 
constant coefficients and for R (1691) were found in f52, (53] as Gauss hypergeometric function, denoted there 
by f = F(X) := F(xu X 2 , X'3> x )i as (for some constants A and B ) 

F( x ) = AF(xi,X 2,X3', x ) + BX 1 ~ X3 F(xi - X,3 + UX 2 - X3 + 1,2 - X3', x )- 

This provides a proof of the statement that f(R) gravity can indeed describe ACDM scenarios without 
the need of an effective cosmological constant. Working with auxiliary connections of the type D, we can 
generalize the constructions to off-diagonal configurations and various classes of modified gravity theories, 
where A, B and Xii X 2 , X3 are appropriate functions of the h coordinates. For instance, reconstruction proce¬ 
dures for Finsler like theories and cosmology models on tangent and Lorentz bundles, and bi-metric/massive 
gravity models are given in 

Having chosen f = F(X) for a modified gravity, we can go further and mimic an off-diagonal configuration 
when f = f(R) is introduced in (12511 and the source Y(x*, £) = A/ 1 f+f/2 x f is considered instead of (1601) and 
(1611) for 9 2 = A^ 1 [$ 2 | Y| + J d( $ 2 ^| Y|], Nevertheless, recovering nonhomogeneous modified cosmological 
models cannot be completed for general re-parameterized dependencies on (x\£) (in particular, only on (). 
This distinguishes explicitly the modified gravity theories of type f(R) from those generated by nonholonomic 
deformations. For certain homogeneity conditions, we can state an equivalence of some classes of gravities 
and cosmological models, or analyze their alternative physical implications. But a complete recovering is 
only possible if all generating and integration functions and the effective sources are correlated with certain 
observable cosmological effects and further approximations and re-definitions in terms of constant parameters 
and functionals depending on a time-like coordinate can be effectively performed. 

In general, a modified gravity theory is not transformed completely into a nonholonomic off-diagonal 
Einstein manifold; an overlap between certain classes of solutions and cosmological and quantum gravity 
models may exist (see the constructions and discussion in mmmmi related to dang ei up. a 
rigorous theoretical analysis of various types of classical and quantum corrected solutions is necessary and 
new experimental data are compulsory in order to conclude that an orthodox paradigm with nonholonomic 
off-diagonal sophisticate (non) vacuum configurations in GR may be enough for elaborating a final, viable 
cosmological model and to perform a variant of the effective covariant anisotropic quantization, what is 
indeed missing, up to now. In a more radical case, we will have to modify substantially the GR theory and 
a number of additional issues may arise on the status of off-diagonal solutions, on methods of quantization 
and recovering formalism, on stability issues and some other, in the search for a matching solution. 

4.2.2 Nonholonomic configurations mimicking phantom and non-phantom matter in f(R) 
gravity 

The anholonomic frame deformation method allows to reconstruct off-diagonal configurations modeling 
f(R) gravity and cosmology with non-phantom or phantom matter in GR. With respect to N-adapted frames 
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in an off-diagonal (modified, or not) gravitational background, the FLRW equations can be written as 

3 k~ 2 H 2 = s p{x k )a~ c + P p(x k )a< xk \ (72) 

where a(x k , Q and H(x k , </) are determined by a solution (1551) . or (15911 . For re-parameterizations or approxi¬ 
mations with (a/, C) — >• C, we assume that the positive functions s p(x k ), P p{x k ) and c(x k ) can be considered. 
The first term on the rhs dominates for small a in the early universe, as in GR with non-phantom matter 
described by an EoS parameter w = —l + c/3 > —1. Introducing q(x k ,C) '■= H 2 (x k ,C) and the respective 
functionals s q := ^ s po-Q C and p q := ^ p pa q, for q = s qe~ cCy + p qe c( =, in R = 3 d^q(Q + 12g(</), we find 


e cC = 


[R ± \Jr 2 -4(144 -9c 2 )]/6(4 + c), 
R/24, 


for c / 4; 
for c = 4. 


(73) 


The non-phantom matter may correspond to the case c = 4 in (1 73 [1 . including radiation with w = 1/3. Eq. (1721) 
transform into a functional equation on Y determined by changing the functional variable R 2 = —576 s q p q Y, 
4Y (1 — Y)-j^j + (3 + T) — 2f = 0. This is again a functional variant (if we consider dependencies on x k ) 
of the generating Gauss’ hypergeometric function, similarly to (1711) . which can be solved in explicit form. 

For the case c / 4 in (I73|) . we come to models with phantom-like dominant components. A similar 
procedure as for deriving Eqs. (22) and (23) in [52 , 53] , results in a functional generalization of the Euler 
equation, namely 


dR 2 


+ AR 


df(R) 

dR 


+ Rf (R) = 0, 


for some coefficients A = —Hq(1 + Hq) and B = (1 + 2Hq)/2, for Hq = 1/3(1 + p hW ). Here we consider, 
for simplicity, homogenous limits and approximations H 2 (t ) = ^ ph p for the phantom EoS fluid-like states, 
phP = phU> phPi with p hW < —1. In both cases, with a trivial or a nontrivial nonholonomically induced 
torsion, there are solutions of the nonholonomic Euler equations above which can be expressed in the form 
?(R) = C' + R m+ + C- R m_ , for some integration constants C± and 2m± = 1 — AY ^/(l — A ) 2 — 4 B. This 
reproduces with respect to N-adapted frames the phantom dark energy cosmology with a behavior of the 
type a(t) = ao(t s — where t s is the so-called Rip time. If the generating functions for the off-diagonal 

cosmological solutions are chosen in a way such that the N-connection coefficients Wi and rq transform to 
zero, the solutions describe universes which end at a Big Rip singularity during t s . Additionally to the former 
result that in the f(R) theory no phantom fluid is needed, we conclude that for off-diagonal configurations 
we can effectively model such locally anisotropic cosmological configurations. 


4.2.3 On nonholonomic constraints and non-conservation of the effective energy-momentum 
tensor 

One can encode and effectively model various types of cosmological solutions for modified gravity theories 
with f(R) and/or f{R,T,R a pT a P) functionals and their nonholonomic deformations. The cosmological 
reconstruction procedures can be elaborated for various types of viable modified gravity which may pass, or 
not, local gravitational tests and explain observational data for accelerating cosmology, dark energy and dark 
matter interactions [23] [25] [26, [52] [53], 32] [33] [33, 35] [46] [37]. Nevertheless, these theories exhibit certain 
specific problems as non-conservation of the energy-momentum tensors for the effective or physical matter 
fields. 

Let us discuss the “non-conservation” issue which is related to the nonholonomic deformations of GR. 
Even in the case V —> D = V + Z, we have the condition D 7 T a ^ ^ 0, which is a typical one for nonholonomic 
(subjected to non-integrable constraints) mechanical or classical field theories. In Lagrange mechanics, for 
instance, the issue of nonholonomic restrictions is solved by introducing additional integration constants. 
In such cases, the conservation laws should be re-considered by taking into account various classes of non- 
dynamical functions. Because the distortion tensor Z is completely defined by the data (g.N), we can 
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compute in unique form the value D 7 T Q( g and relate this to the fact that, in general, R a ^ is not sym¬ 
metric. This is a consequence of a nonholonomically induced torsion. It is convenient to work with such 
nonholonomic variables (g, N.D) in order to apply the anholonomic frame deformation method and decou¬ 
ple certain modified gravitational equations and generate off-diagonal solutions of (1231) . After this general, 
integral variety of solutions has been found, one can re-define the generating functions and sources in order 
to generate LC-configurations, as was shown in section 13. 1.3 1 This way, the problem of “non-conservation” of 
effective and physical T a g can be solved by encoding into generic off-diagonal configurations with effective 
conservation laws which are similar to GR. 

In explicit form, we explain how the “non-conservation” problem can be solved for off-diagonal solutions 
with one Killing symmetry in the framework of /(/?,, T ) theories generalizing certain constructions from 
j54| . Following a similar procedure as in Sect. II of that work, but using the operator D instead of V, for 
f = f(R, T), and considering an N-adapted parametrization of the effective source T = const , we prove 
that 


(1 + ^)D“T af} = Jg a/ 3 D“T - (T q/3 + 0 q/ 3 )D“ ln( 2 f) - D“0 a/3 . 


(74) 


In these equations the values 2 f := df/dT and Q a p = — 2T 0( g — pg a /3 are used, with an energy-momentum 
tensor (I26p for nonholonomic flows of a perfect fluid. In general, D Q T aj g ^ 0 even for nonholonomic defor¬ 
mations of GR. Nevertheless, we can consider a subclass of off-diagonal configurations in f(R, T) gravity 
when T (1321) — > A (1331) and 4 > 2 = A _1 [<I> 2 | Y| + f d( <h 2 3^| Y|] (16T1) result in f — > f = R and effective 
R°/ = Ad 0 / which admit LC-solutions with zero torsion. For such nonholonomic distributions with D —>• V, 
D“T q/3 —» VA = 0 and all terms on the lhs of (1741) vanish, because they are nonholonomically equivalent 
to functionals o f the effective cosmological constant A. Such conditions are satisfied in correspondingly N- 
adapted frames and for canonical d-connections. The equations (j74[) generalize to nonholonomic forms the 
similar ones derived for the Levi-Civita connection V (see Eq. (10) in Ref. [54]). 

The already mentioned problem of “non-conservation” becomes worse for general / (R,T , R a pT a P) the¬ 
ories. Even if there are certain special cases where it can be solved [241125 , 2611521 (S3], it is the case that 
no solution can be found in general form. Surprisingly, the anholonomic frame deformation method also 
suggests a procedure for selecting off-diagonal configurations which can mimic modified gravity theories in a 
self-consistent way, admitting effective torsions completely determined by observational data (f, Y, g, N,D), 
and/or by constraints to LC-configurations. By redefining the generating functions in the form Y (1321) —>• 
A (I33[) . we get the possibility to consider effective sources and off-diagonal Einstein spaces for which the 
divergence and non-conservation problem become similar to those in GR or (for more geometrically complex 
configurations) to those in viable f(R) models. Finally, we emphasize the fact that in this way we do not 
find a general solution for all f(R, ...)-modified theories but only for those models which admit an encoding 
in an effective GR system and a general decoupling via the nonholonomic deformations. 


4.3 Nonholonomic constraints and matter instability 

There is another serious problem in modified gravities which is related to possible matter instabilities 
related to modifications of the gravitational actions. Even tiny modifications of GR may make the new 
model to posses unstable interior solutions (see, e.g., [55 . 56]). It was demonstrated however that there are 
viable f(R) theories (with appropriated choices of the functional) where such instabilities may not occur 
dU2U3uiiEi[6i[7U8UHi uni [m na na nn etj. In this section, we apply the anholonomic frame deformation 
method to more general f(R,T,R a pT a P) theories. The corresponding field equations are very difficult to 
solve even in a linear approximation [241 [25] [26] . if we work in coordinate frames and general functionals. 
In the nonholonomic variable formalism, the gravitational field equations in modified gravity theories posses 
the decoupling property exhibited above, which allows to encode f(R, ..^-modifications into off-diagonal 
nonholonomic configurations for the effective Einstein manifolds. 
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For a stability analysis, the trace equations where (1311) are multiplied by are to be considered, namely 


- 2? + (R + 3D At D M ) *f + (T + 0) 2 f + XD^D^T + D #1 D I/ T'“' + S) 3 f = n 2 T, 


(75) 


where 1 f := <9f/<9R, 2 f := di/dT and 3 f := <9f/<9P, when P = R a ^T a/3 . Let us envisage a trace 
configuration in the interior of a celestial body, when T = To and — 2f + Ro ( 2 f) = k 2 To- Imposing 
nonholonomic constraints, we parameterize a LC-configuration in GR and model an interior solution in the 
presence of some gravitational objects (for instance, the Sun or the Earth). The /-modifications (in general 
with strong coupling for the curvature and the energy-momentum tensor) may result in a worsening of the 
stability problems and may prevent To to be a solution of any suitable background equation. It is difficult to 
find solutions of (1751) even for very much simplified cases in the nonlinear situation if we work in coordinate 
frames for the connection D = V. 

A rigorous study of the problem of matter instability for f(R) and more generally f(R,T,R a pT a P) 
gravities, for certain illustrative cases when 1 f = R, and for restrictive conditions where there is a qualitative 
description via additional functionals on T and P shows that the appearance of a large instability can actually 
be avoided. Using the anholonomic frame deformation method, we can consider modified gravity theories 
with /-modifications which are effectively described by R°jg = Ad^ when the modifications are encoded into 
polarization functions and N-coefficients. For models generated by 

f(R, T, R q/ 3 T^) = ?i(R) + F(P) + G(T), 


we take a constant interior solution with To = const and Po = const , and denote by f|' s ' 1 := <9 s fi/clR s and 
F(s) := 

d s F/dP s for s = 1, 2,... We can repeat, with respect to the N-frames (1121) and (1131) . the computations 
presented in detail for Eqs. (45)-(48) in |2411251126| (see also references therein), and prove that Eqs. (1751) for 
linear pertubations can be written in the form 


.Po fi 


(i) 


[D^Du + 2^2—D, t l\ + 2^ + 4-. „ 

L ^ T 0 M T 0 T 0 F(2)- 


ip = [, 2 f > 


(i) 


Po FW 


T 0 F( 2 ) T 0 F(2) 


X-T n 1 SR. 


The values of type 4P and dR are considered for a small perturbations in the curvature where R = Ro + dR 
and P = Po + dP- No instability appears if dP =<7R = 0 which is a particular solution of the above equation. 
We can in fact model a damped oscillator with additional nonholonomic constraints if Ho+2Pof{ 1 Vf ( 2) > To, 
which allows to avoid large instabilities in the interior of a spherical body. For some specific functionals f(R) 
and appropriate G(T), the same behavior as in GR results (with mass stability in the cosmological context), 
although there are possible large perturbations SR. and SP remaining. The ideas how to circumvent the mass 
instability problem for holonomic configurations has been studied in | [58| [59 , 60] [6TI [62 ,1 63] |6~4] . Redefining 
the generating functions and sources in a /-modified model into an effective Einsteinian theory, with S [R,A], 
one can consider a nonholonomically deformed Hilbert-Einstein action with f —» f = R. In such cases, 
e)R = <7R = 0 and instabilities are not produced, indeed, if we impose the zero torsion conditions (see (1241) 1. 
we get back to the GR theory. Even if Eq. (1751) involves not only perturbations of the Ricci scalar R but also 
of the Ricci d-tensor Rq.^ (through <5P), via nonholonomic transforms to effective R°g = A d a p, the stability 
of the system is obtained via off-diagonal interactions and the nonholonomic constraints used for an effective 
modeling of a subclass of f-theories to certain nonholonomic deformations of the Einstein equations with 
effective cosmological constant A. This is indeed a remarkable result. 


5 Effective field theory for off—diagonal cosmological configurations 

Both in particle and condensed matter physics, effective field theory (EFT) methods have proven so far to 
provide a quick and economic way in order to connect experimental data and phenomenological results with 
certain fundamental theories (see the reviews |651166] ). Recently methods of that sort have been applied in 
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cosmology, in particular to inflation [1571 [68] 69], late-time acceleration (SUEDES and dark energy physics 
[76. iTl] (details and references can be found in |75l 76]). The goal of this section is to construct an EFT 
describing perturbations both over diagonal and off-diagonal cosmological background solutions in modified 
gravities, in cases where an effective Einsteinian manifold can be associated, and when the matter sector 
obeys the weak equivalence principle and all modifications of gravity and the matter fields can be encoded 
into an effective cosmological constant. 

5.1 Off-diagonal background evolution and ACDM 

We shall consider here configurations with redefined generating functions and sources and where the third 
effective action in (125[) is taken to be of the form 

S= 9 S+ m s = J { _ i_ fi(t)R + A(t)- c (t )( i g 44 ]^id 4 u+ m S[g a/3 ]. (76) 

This is related to a background FLRW solution constrained, on its turn, by observation to be in close 
agreement with ACDM^fl In this formula, the value dg 44 is the perturbation of the time—time like component 
of the d-metric, and the Ricci scalar R is defined by a canonical d-connection D which can be obtained by a 
finite chain of redefinitions, resulting in [g, N,D, T] —7 [g, N, D, T a g, A(i)] where A(f) takes a constant value 
A/27c 2 for the /i-equations while it is a function of t for the u-equations. The energy-momentum tensor 
T a /3 determined by m S[g a ^] does not satisfy, in general, the condition D Q T“^ ^ 0 , but V Q T Q ^ —> 0 for 
D„ —>• V Q . For simplicity, effective matter is treated as a perfect fluid, which in the N-adapted model is 
described by a time dependent average energy density ~p(t) and pressure p(t). For LC-configurations one has 
the usual continuity equation 

-f + 3 H{p+p) = 0, 

where H is determined by a scaling factor a(t) for a generic off-diagonal solution of type (1551) . or (l59l) . We 
note that in this section p 0 is the derivative with respect to the physical time (the formalism works also in 
conformal time). The function a{t) can be considered as the limiting result for e —>• 0 of an integration in a 
metric of type Q (where nonlinear interactions are encoded by f-modifications), as we explained before for 
that formula. 

5.1.1 Off—diagonal background evolution 

The background evolution is in general off-diagonal (values e ~ 10 -20 do not contradict present exper¬ 
imental data [5l]). For simplicity, we can consider a diagonal background with a FLRW metric with zero 
spatial curvature. We can chose a(t),~p(t) and p(t) to be close to ACDM evolution if e —>• 0. As in standard 
(diagonal) cosmology we can use the Friedmann equations to eliminate the functions A (t) and c(t) (which 
can again be considered as a redefinition of the generating functions in our approach) but keep Q(t) as a 
free function, similarly to |I3 [76] and [77]. For diagonal configurations, such theory can be formulated in 
the Jordan frame and dealt with a non-minimal coupling between an effective scalar field and metric. In the 
Einstein frame, we get a coupling of matter to the effective scalar field. Nevertheless, the scaling factor a(t) 
has some “memory” of the genuine nonholonomic redefinition of the integration functions and corresponding 
contributions of the modified gravity theory. 

Varying (176[) with respect to N-adapted frames, we get an effective Friedman equation, which allows to 
express 


-2 c(t) = p + p + (2nH°+ n 00 /2-HQ*)/*; 2 
-A (t) = p + (3 nH 2 + 2QH + + R9°) /7c 2 

15 In effective field theories the mass scale to 2 = k -2 can be different from the Plank mass of GR. It is used to render ff 
dimensionless. 


25 



for any given data (fl, o, /0,p)El These equations describe how the background energy and pressure of the 
DE component evolve over cosmological history, corresponding to the evolution of the N-coefhcients Wi(t). 


5.1.2 N—adapted perturbations 

Assuming that for D —>• V the weak equivalence principle is true, one can always introduce a conformal 
(Jordan) frame, when the matter fields couples only to the d-metric g a p and not to the scalar field. 

Starting from an N-adapted form for (1 76 [1 . a procedure similar to that for the construction of effective 
field theories (see Eq. (2.1) of |751176] 1 leads to the general effectively unitary gauge action 


5 = 


I + A(t) - c(f)<5g 44 + ^W 4 ) 2 + ^M(hg 44 ) 3 + ... 


(77) 


M Hg 44 JK« Q -4^K^_ Ml{t) 


5K a *5KP a + ... + 4 A(i)(hR) 2 + 2 A(i)hR a/3 hR a/3 


2 ° 2 2 ? 

+ 1 7 (t)C^ (jA C^ A + 2 7 (t)e^ A C^ aP C aXa p + ... 

+ 1 m 2 (t)n M n !/ (e /U g 44 )(e iy g 44 ) + 2 m 2 (f)( g^ + n^n")(e /i g 44 )(e !/ g 44 ) + ...] + m S[g af} ]. 


The value C^ LUrT \ is the Weyl tensor determined by D, while rh* is the vector normal to the surfaces of 
constant time. Each term in (|77|) can have a time-dependent coefficient because the background solutions is, 
in general, off-diagonal and breaks time translation symmetry. The matter action m S[g Q ^] can be arbitrary, 
with sources to be re-defined by integration functions. We will fix below the quantities <5g 44 , <5R,hR a rt and 
hK a p as perturbations, the dots corresponding to various terms which we do not specify of quadratic and 
higher-order perturbations. N-elongated partial derivatives (11211 are used instead of d^. 

We involve and additional 3+1 splitting (the variant 2+2 is convenient for constructing off-diagonal, 
exact solutions), where the Ricci d-tensor of D decomposes as 

R Q/3 n Q n^ = K q/ 3 K^ - K/K/ + D Q (rED 7 n Q ) - D Q (n“D 7 rE). 

As a result of the N-adapted construction, the value n Q JK Q , ( g, where <5K a p = K a/ 3 — °K Q , / 3 := K a p + 
3H(g fJ,u + n^rE) vanishes. 


5.1.3 Effective field theory in terms of Stiickelberg d—fields 

In effective held theories and cosmological models, the Stiikelberg technique (SZl Eg ES Eg ESI is used 
when explicit functions of time are modified according to 

ip(t)ip(t + ir(u a )) (78) 

and then Taylor-expanded in ir(u a ). Such procedure is applied to the action (1771) when the ansatz for off- 
diagonal solutions (1551) , or (159(1 , is reparameterized as a 3 + 1 form, 

ds 2 = g a/3 du a du 13 = a 2 (gjj + gij)dx-dx 3 - — dt 2 , synchronous gauge ; (79) 

= a 2 (gij + <fy) deeded- — (hi) 2 , N-adapted synchronous gauge ; 
ds 2 = a 2 ( 1 — 2cj))gijdx L dx 3 - — (1 + 2iJj)dt 2 , Newtonian gauge ; (80) 

= a 2 (l — 2<j))gijde L dei- — (1 + 2'0)(hf) 2 , N-adapted Newtonian gauge , 

16 This can be rewritten in a more conventional form in terms of the dark energy density, pde, and pressure, pde, (cf. (|4])1. 
when 

3 QH 2 = k 2 (pdb + p) and 2f IH° = —k 2 (pde + Pde + p + p), 
for 2c(f) = pde + Pde + - f+°)/7c 2 and A(t) = pde - (tffi 0 + Q. 0<> )/k 2 . 
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where de- = (dx l , 5e 3 = dy 3 + hidx l ) and St = dt + Widx\ see (USD. The basic premise for this is that we can 
perform coordinate transformations t = t + ir(u a ),x- = x -, were the convention for indices is i,j ,... = 1, 2,3 
and the coordinates u a = ( x l ,y 3 ,t ) —>• u a = (x-,t). We consider to be a time-independent background 
metric but a{t) and [N °‘. g t , h a \ are certain nonlinear functionals determined by a cosmological off-diagonal 
solution in the modified theory. The metric (1791) is written in a form which allows to describe perturbations 
in a synchronous gauge. This parametrization can be obtained for any d-metric m which with respect 
to coordinate frames is rewritten in the form (1751) with coefficients (1 16 [) . In the synchronous gauge, the 
coordinate transforms u a —>• u a are chosen in a for to be satisfied the conditions Sg 41 = g 4 - = 0 (we note 
that here we do not use boldface symbols because, in general, such conditions can be imposed for not 
N-adapted frames). It will be convenient to discuss some phenomena in the Newtonian gauge (1KU1) . For 
small e-deformations (see ©), off-diagonal extensions of cosmological metrics can be treated as effective 
fluctuations which may be nonholonomically constrained, or not, to be parameterized in N-adapted form. In 
both cases, we can use the formulas derived in [73176] but keeping in mind that we work with a d-connection 
D which only at the end will be additionally constrained, when D —> V0 

We can decompose 

+ (DjD i - ^D 2 ) 7 ? and Qj = ^g^ + (Dpj - ^D 2 )? 7 , (81) 

where tildes refer to quantities associated with the spatial metric, and express the metric determinant *J\g\ = 

a 3 \J TM(1 + ?/2) and ^/jg[ = a 3 y / [§|(l + <f/2). It is convenient to use a definition in mommentum space 
k a = (k ,k), where 

rf(t ,t) = —k~ 2 [g(lz , t) + 6 rj(t,t)] (82) 

(see details on similar conventions in [75; 176]). where the effective held constructions are also performed in 
the massive case. 

5.2 Linking off-diagonal perturbations to observations 

The goal is here to identify and study the full set of perturbations which result in off-diagonal deforma¬ 
tions of the FRW background. Linearized equations of motion for certain effective scalar perturbations will 
be considered, which will allow to determine: 1) the speed of sound; 2) Poisson’s equation; 3) the anisotropic 
stress; 4) the effective Newton constant; and 5) Caldwell’s parameter [78] . Reparametrizations of the metric 
in the respective forms (1791) and (1801) allow for this association with the standard phenomenological functions 
and parameters appearing in literature, connecting in this way observational features with the off-diagonal 
cosmological solutions, for different modified gravities and effective Einstein spaces. In what follows, and 
as a hint to explicitly show these possibilities, we will consider in brief 1), 4) and 5) above, as obtained for 
certain classes of cosmological models. 

5.2.1 Effective speed of sound in off—diagonal media 

The data and fj defined by Eqs. (1781) . (1811) and (182p are used for writing the so-called 7 r-equation of 
motion along with the time-time and space-space trace components of the Einstein equation (in our case, 
modified in terms of V —> D) via the kinetic matrix 7 x,y (namely, the coefficient of Y in the X equation 
of motion). In its turn, this allows to compute the speed of sound of scalar perturbations in the sub¬ 
horizon limit using the synchronous gauge for this calculation. The differential operators are transformed 

17 For simplicity, we do not dub, in an N-adapted form, the proofs from those works but use directly the synchronous gauge 
representation which is more convenient for studying both perturbations and e-deformations all included in a term Q, (e. x L , l); 
we omit the tilde on spacetime coordinates when this does not result in ambiguities, and use in brief the term perturbations both 
for e-deformations and for fluctuations; in a more general context, perturbations with respect to a fixed N-adapted background 
can be considered. 
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into matrix components by changing z* —>• —ico. This allows to compute the coefficients of CZZD, following the 
same procedure as in Table 1 and Appendix D (see also Eqs. (4.1)-(4.5)) of [73, 176| . where the determinant 
of the kinetic matrix is set to be zero. Such equations are parameterized in the form 


/ 77r,7r 7 nw Tn-,17 \ / 7 T \ 

7 tt,n 7 u,c Mv l = 0 

V 1ss,ir 7 ss,<; Isstf / \ V J 


The resulting expression yields, in general, a number of nonlinear dispersion relations. For instance, the 
presence of the operator M^{t) results in the dispersion formula for the speed of sound, c, 


Cs 2 = 1 + 2 M 2 4 


c + 


3(17 


o\2 


47t 2 17 


-1 


The value c is linked to various stability issues, as analyzed in [721 EH, where it was concluded that ghost-free 
conditions are satisfied if 


2 


3(17*) 2 1 

c+ lFiT 


PDE + PDE + «' 


HiY 


, 3 (17°) 2 

2 n 



> 0 . 


For certain off-diagonal configurations, the effective EoS of DE can become of phantom type, wde < “’ll 
even if the nonholonomic deformations/perturbations do not seem to host one (it may depend on the choice 
of generating function 17(t)). 


5.2.2 The effective gravitational constant and Caldwell’s parameter 

In the Newtonian limit, we look at off-diagonal deformations and modified theories of gravity that share 
the horizon k 2 /a? H 2 . Such a a horizon is indeed shared by different types of solutions in different models, 
and may show up, for instance, in e-deformed nonholonomic backgrounds. Similar considerations as for 
subsection 4.4 in [231251 result in an expression for the effective Newton constant 47T e ^G := k/2 (note 

here 7 r is the mathematical constant, and not the function ir(u a ) from (178p ) 


e ff K 2 = K 2 n- 1 


1 + (17*) 2 /4 17(c7t 2 + 


2 3 (fl°) 2 2 ^ 2 t, r2\ 

■^ 2 1 v ; +a 2 —M 2 ^ 


4 n 


k 2 ' 


This expression becomes more complicated if more terms, in addition to M 2 , are included. There are 
modifications of G even if 17* = 0. We note here the explicit dependence on a 2 , which can be considered 
as the generating function for off-diagonal solutions. 

Let us now consider the parameter c w considered by Caldwell et al. m as the cosmological analog of 
Eddington’s PPN parameter 7 , when c w ~ 1 — 7 (we use here a notation which differs from those in other 
works on the subject, in order not to interfere with zu from Eq. ([3]), etc.)@ Using background terms, one 
obtains 


VJ = 




o\2 


217 


cTc 2 + ( °' 


o\2 


217 


+a 2 i Afz 


1 -1 


Such formula becomes more complicated if further terms are included. In both cases of this and the previous 
formulas, only the terms M 3 and M 2 = (5g 44 (fR ( '^ (for perturbations of the spacial scalar curvature R® 
of D) contribute, if 17° = 0. Both n 2 (a) and c w(a) are functionals of the generating functions. In GR, 

e fftz 2 = 1/mp and c w = 0. Nevertheless, these values are nonzero if the configurations are off-diagonal with 
an effective cosmological constant. A kind of nonlinear classical polarization of the gravitational constant is 
there obtained at certain scales. 

ls There are bounds of the type c zu = 0.02 ±0.07(2 a), for scales of about 10 Kpc, and =0.03 ± 0.10(2c), for hundreds of Kpc, 
but as of now there are not known limits at the Mpc scale yet. 
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5.2.3 Confronting an off diagonal cosmology which encodes modified gravity with actual ob¬ 
servational results 

In Sect. 4.6 of [751 76] (see, in particular, Tables 4 and 5 there), the contributions from various operators 
in a generalized model with an effective field theory action were calculated. We will here combine this 
knowledge on the dark energy models (in our case with off-diagonal interactions modeling the ACDM, /(R), 
and Horava-Lifshitz |79j theories) for the nontrivial N-adapted coefficients in (17711 : 


Operator 

n 

A 

c 

M 2 

2 m 2 

MGTs 

R 


<5g 44 

m a a ? 

a _2 g^(e i g 44 )(e j g 44 ) 

ACDM 

1 

Y 

0 

— 

— 

/(R) 

T 

Y 

0 

- 

- 

Horava-Lifshitz 

1 

Y 

0 

Y 

Y 


In the above matrix, the symbols for the respective operators mean that T is necessary, the — is not included, 
and 1,0 is for unity or exactly vanishing, respectively. 

We conclude that off-diagonal modeling of modified gravities can be performed as small perturbations 
of ACDM, when dark energy exists in the form of a cosmological constant. The off-diagonal redefinitions 
of generating functions result in a modification of a(t) which is different from that in diagonal cosmological 
models. Nevertheless, even very small values of e-off-diagonal deformations may introduce a certain speed 
of sound and clumping of DE (this is different from diagonal configurations in when such effects are zero 
and the effective Newtonian constant is just that for GR). Nevertheless, the models can be characterized by 
different physical values if the off-diagonal modeling is performed for different modified gravities, indeed: 


Operator 

n 

M 2 2 

2 m 2 

MGTs 

R 

m a Q ) 2 

a' 2 g^(e i g 44 )(e i g 44 

Speed of sound 

1 

Y + k 4 

* 

Effective Newtonian constant 

Y 

Y + k 2 

Y 

Caldwell’s parameter c vo 

Y 

Y 

— 


Here additional labels have been used, for instance, Y + fc 4 is for a new scale dependence, and * is set for an 
operator behaving unusually. 

Finally, we note that the effective field theory can be generalized in a form which allows to describe 
dark energy and modified gravities in the late universe when off-diagonal nonlinear parametric gravitational 
interactions encode contributions of such modified theories. The action for perturbations (1771) of the effective 
models (1761) includes positive small e-deformations of the ACDM models with nonholonomically constrained 
fluctuations which contribute to the background cosmological evolution. A systematic investigation of the 
physical effects of N-adapted operators was undertaken which emphasizes the importance of off-diagonal 
effects when modeling modified gravities with the aim to constrain the possible nature of DE. 

6 Short summary of new geometric ideas and methods for the construc¬ 
tion of off—diagonal solutions 

In the modern era of precision cosmology, a series of evidences have been found pointing towards devia¬ 
tions from the standard big bang cosmology. The universe may be slightly non-homogeneous and anisotropic 
at very large scales and the accelerating expansion phases are determined by dark energy and dark matter 
effects. A number of alternative modifications of gravity have been proposed already, with the aim to elabo¬ 
rate realistic cosmological models. Some of the most intensively exploited are the so-called f(R, ...)-modified 
gravity theories HlElElllUSlEUTlISUSlIinilllUISUISUIll- The gravitational field equations in GR and the 
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modified gravities consist of very involved systems of nonlinear PDE. A rigorous study of the cosmological 
evolution in the frames of the different theories requires new analytic, geometric and numerical methods for 
constructing exact and approximate solutions, reconstructing procedures and effective field models. In sce¬ 
narios closely related to the standard ACDM universe, the ansatz for the alternative metric is usually taken 
of FLRW diagonal type and the interactions are modeled by effective matter, exotic fluids and modifications 
of GR. The advantage of the anholonomic frame deformation method, AFDM, [2411251 [26] is that it provides 
a geometric formalism for constructing exact generic off-diagonal solutions encoding nonlinear parametric 
effects which mimic scenarios of anisotropic cosmology via generating and integration functions. Such nonlin¬ 
ear gravitational interactions, and the possible cosmological implications thereof, are not considered if we fix 
from the very beginning the diagonal ansatz, which results in a system of ordinary differential equations. For 
generic off-diagonal configurations, we can always derive an effective field theory and confront cosmological 
theories with existing observational data as we showed in Sect. [5] 

For the equivalent modeling of exact solutions and cosmological scenarios in different classes of modified 
and GR theories, we consider three different parametrization of the action for gravitational and matter fields 
(1251) . with respective Lagrange functionals 

C = g jr[f(R,T,R a 0 T^)] + m £[g a p,V, m <p\ (83) 

= 9 L[f(R, T, R Qj gT“^)] + m L[g Q ,^,D, m (p\ (84) 

= g t+ m L = R+A. (85) 

In these formulas, are some (effective) matter fields, which can be approximated by the components 
of perfect (pressureless) fluids, for instance, with an energy-momentum tensor of e thtype (126[) . and the 
linear connections V, D and D are related via distortion relations of the type (1191) . which are completely 
determined by the metric field g a p ~ g a/3 ~ g a/ 3 up to frame transformations. The functionals 9 C[f{...) ] 
and m C[..] determine the corresponding model of /-modified theory of gravity. The Lagrange densities are 
written in different geometric variables because ofthis allows us to find exact solutions and model physical 
effects by means the same solutions but in different theories of gravity. 

6.1 Decoupling of the generalized Einstein equations and modified cosmological solu¬ 
tions 

We can decouple the gravitational and matter field equations (l3ll) in a modified theory of gravity, derived 
from a Lagrangian density (I83[) . and construct generic off-diagonal solutions depending (in general) on all 
spacetime variables, if we work with geometric data (g, D, N). The generalized Einstein equations are written 
in the N-adapted form R Q( g = Y aj g with an effective source X Q;j g = diag[ h T, ”Y] determined by the matter 
fields and the nonholonomic and /-deformations. For a very general off-diagonal ansatz for metrics with one 
Killing, or non-Killing, symmetries, the gravitational field equations transform into a system of nonlinear 
PDEs (I40[) (14411 . which can be solved in their general form for nonhomogeneous and locally anisotropic 
cosmological configurations (1551) . In the coordinate frame u— = (x k ,y 3 ,y 4 = t) (where t is the time like 
coordinate and k = 1,2), such metrics ds 2 = gapdu—dvP- are of the type © , or ([29]) . and can be written in 
a form involving a generalized scaling factor a 2 (x k ,t), 


a 2 r ]i + u} 2 [(ni) 2 a 2 h 3 - (uq) 2 ] 

oj 2 [nin2a 2 h^ — W 1 W 2 ] 

uj 2 nia 2 hs uj 2 w\ 

w 2 [nin 2 a 2 /i 3 — uquq] 

a 2 r ]2 + {n 2 ) 2 a 2 hs - (w 2 ) 2 

u 2 n2a 2 hs u 2 W2 

uj 2 nia 2 h 3 

W 2 7l2a 2 /l3 

u} 2 a 2 h,3 0 

U 2 Wl 

U 2 W2 

0 -uj 2 


( 86 ) 


By fixing the data for the generating functions: <F(x fc , t), 4 > 0 = c/<f> ^ 0; u 2 (x k ,y 3 ,t) as a solution of 
diOJ — (c/^/d^uA = = a 2 (x k ,t)gi(x k ,t) = a 2 (x k ,t)r] 2 {x k ,t) as a solution of ip** + ip" = 2 h T , 

for the nontrivial sources h T(x k ), v T(x k ,t) and for the effective cosmological constant A / 0 , we can the 
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express the coefficients of this metric in a certain general form. One further obtains [see the related formulas 

(1501) . (15TT) . m aR d (1551) ] h-s = /i- 3 /a 2 |h 4 |, with 


/i 3 (.x fc ,t) = <I> 2 /4A, h ,4 = a 2 (x k ,t) = 


(d> 


2\o 


<b 2 | ”T| + / dt $ 2 | ^T|° 


Wi (x k ,t) = $$[8, w T]/$°[$, *T], n k = ink + 2 n k J dt h A /{^\h^\f 


-l 


where in k (x l ), 2 n k (x l ) are integration functions. The solutions (I86[) include, in general, a nonholonomically 
induced torsion (l24l) . which is important for constructing exact solutions in modified theories with nontrivial 
torsion fields. 

We can always put additional zero-torsion constraints (1391) and construct Levi-Civita configurations. 
Such solutions are extracted by choosing 2 n k = 0 and 4 n k = d k n with a function n = n(x k ) and for a 
subclass of generating functions for which (c),;^) 0 = (9*$°, see the details for Eq. (1561) . One 

finds Wi = di&/& = d{A for a nontrivial function A(x k ,y A ) taken to be a solution of a first order PDE 
with effective a sources depending functionally on 4?. The u-components of geometric and physical object 
are generated by couples of data (<f>, ^Y) and (4>,A), or related by formulas of the type (l49l) . or (for zero 
torsion), of the type (1611) . for a fixed value of the effective cosmological constant A = A, with 


A6 2 = $ 2 | W T| + 


dt $ 2 | ?, Tr 


(87) 


For the data (4>[4>, U Y],A), a metric (l86l) is equivalent to the d-metric (1551)1^1 

The fact that the AFDM setting allows us to integrate in this general, off-diagonal (with Killing or 
non-Killing symmetries), both the gravitational field equations of GR and of modified theories, is certainly 
a very important result in mathematical relativity. This result has also fundamental physical implications 
for modern standard and modified gravity theories, particle physics and cosmology. The first one is that 
off-diagonal solutions of (generalized) Einstein equations, depending generically on three or four spacetime 
coordinates, can be generated by general classes of generating and integration functions. This reflects a 
specific property of nonlinear and nonholonomic off-diagonal gravitational interactions, where re-definitions 
of generating functions (for instance, of type (I87[) j of modified gravity theories with sources y a p [see (132 p 
and (1331) ]. allow us to describe a large class of these modified gravity interactions as effective Einstein spaces. 
The effective equations R°g = A<5°g are derived for the gravitational Lagrangian (185 P with a correspondingly 
nonholonomically deformed linear connection D. We have found that we can mimic certain classes of solu¬ 
tions of modified theories as off-diagonal configurations in ordinary GR, and inversely, but using different 
nonholonomic variables and deformed geometric/physical objects. Such constructions cannot be realized at 
all if we chose from the very beginning the ubiquitous diagonal ansatz for the metric and consider holonomic 
configurations. It is not clear yet how to prove the stability of the solution for /-modified theories, with the 
exception of some vary special cases. Nevertheless, for a very general class of such nonlinear nonholonomic 
systems re-defined as f-theories in terms of R°g and R, the stability can be proven indeed, as for the Einstein 
spacetime manifolds (see m • 

Another important physical implication is that, for correspondingly fixed data, for generating and inte¬ 
gration functions solutions of the type g a p(x k ,t) (I86J) determine new classes of nonhomogeneous cosmological 
metrics in GR as off-diagonal deformations of the FLRW cosmology, in the limit g a g —> g a p(t, a(t), hs(t), &(t), 
r}i(t)). We can model accelerating cosmology and dark energy and dark matter effects via nonlinear off- 
diagonal interactions and nonholonomic constraints induced by corresponding transforms of type (1871) and 
modified re-scaling factor a(t) as in d-metrics 055(1 . or Q59[) . The possibility to mimic modified gravities as 
analogous models in GR via nonlinear transformations of generating functions and sources (187[) reflects a 

19 We emphasize that such locally anisotropic cosmological solutions are generically off-diagonal (because, in general, the 
anholonomy coefficients W t b a , see m are non-zero). Some of the six independent coefficients of the metric depend on all 
spacetime coordinates. If we fix u > 2 = 1, we generate solutions with Killing symmetry on 83 . 
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fundamental property of the class of off-diagonal solutions of the gravitational field equations. This prop¬ 
erty holds true for solutions with one Killing symmetry and/or non-Killing symmetries when metrics are 
generically off-diagonal and depend (in 4 dimensions) on three/four spacetime coordinates. It reflects a 
specific nonlinear dynamics of the gravitational and matter fields when off-diagonal interactions are taken 
into consideration for maximally possible six degrees of freedom of the metric and with certain classes of 
nonholonomic constraints imposed. Working only with diagonalizable, two Killing symmetries and stationary 
configurations, such nonlinear gravitational physics and cosmologycal effects cannot be encountered. A rad¬ 
ically “orthodox” interpretation of this class of nonlinear and nonholonomic configurations and cosmological 
evolution scenarios is that they may explain the bulk of accelerating cosmology data and related dark energy 
and dark matter effects. It might be the case that, in order to understand the observed Universe it is not 
enough to modify GR in a simplistic way but, rather, to bring into consideration, within the standard GR 
theory, of a richer class of off-diagonal solutions and then take convenient limits leading to effective theories, 
in the end. This issue still need rigorous theoretical and observational consideration. 


6.2 Small off— diagonal /—deformations and effective FLRW like cosmologies 


We can consider the subclass of generating functions, effective sources and cosmological constant 

U T],A) as in (1H71) . where the cosmological evolution of spacetime regions is approximated by off- 
diagonal deformations with polarization functions rj a ~ 1 + ex a (x k , t) and the N-coefficients rii(x k ) and 
Wi(x k , t ) are proportional to a small parameter e when 0 < £ <C 1. This is motivated by the fact that although 
possible anisotropic cosmological effects are very small, the modifications of the scale factor a(t ) —» a(x k ,t) 
[for a FLRW metric (12711 with for g\ = 92 = 93 = a 2 , <74 = —1] can he substantial for some intervals 
of time when the generating functions are of type (18711 . Prescribing any a(x k ,t ) and the solution e^ x ^ 
compatible with the observations, and fixing, for simplicity, c S 2 = 1, as for the d-metric (1591) . when the 
polarization functions can be approximated as 1 + e\i = a~ 2 e^, arbitrary X3(x k U) but x .4 = k and the 
function /13 = 773 /a 2 , we obtain, up to £ 2 , 


9afi_ — 


a 2 (l + ex 1) + e 2 [(ni) 2 a 2 - (uq) 2 ] 

£ 2 [niri 2 a 2 — W1W2} 

£n\a 2 

£W 1 


£ 2 [nin 2 a 2 — W1W2] 

a 2 (l + ex 2 ) + £ 2 [(n 2 ) 2 a 2 h 3 - (rc 2 ) 2 ] 

£n 2 a 2 

£W 2 

• ( 88 ) 

£n\a 2 

£n 2 a 2 

a 2 

0 

£W\ 

£W 2 

0 

-1 



This class of locally anisotropic metrics is of type (18611 , with Killing symmetry on 83 and small off-diagonal 
deformations on an anisotropy parameter £ which has to be fixed by experimental data. We can consider 
the limit e —> 0 in (1881) but even in such cases we have an anisotropic scaling factor a(x k ,t), or a(t), which 
is different from a the standard a(t ) FLRW one. This is a consequence of the nonlinear off-diagonal and 
nonholonomic gravitational interactions, with generating functions and possible modified gravity sources 
related by transforms (1871) . 

Having constructed a class of LC configurations (188[1 . we can extract certain subclasses of cosmological 
evolution scenarios with generating and integration functions when a(x k ,t ) —>• a(t ) 7 ^ a(t),Wi —>• Wi(t),rii —>• 
const., etc. In this way we reproduce an effective FLRW like cosmology when the scaling factor a(t ) is defined 
not by exotic dark matter and dark energy interactions but by certain off-diagonal gravitational interactions 
which mimic contributions of the modified gravity type. In a more general context, metrics of type (I88[) 
may encode certain nonholonomic torsion configurations if the LC-constraints (139[) are not imposed. Such 
solutions also contain a small parameter e but a scaling factor a(x k ,t ) —» a(t ) 7 ^ a(t ) is generated by data 
(<b[<b, U Y], A) related by formula d49l) instead of (1571) . 


7 Concluding remarks and discussion 

We have proven in this paper that a wide class of f(R ,...) modified gravity theories can be encoded 
into effective off-diagonal Einstein spaces if nonholonomic deformations and constraints are considered for 
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the nonlinear dynamics of gravity and matter fields. A special attention has been paid to a new version of 
modified gravity theory which includes strong coupling of the fields [23l:.25,:26]. Such modified gravity theories 
have physical motivations from the covariant Horava-Lifshitz like gravity models, with dynamical breaking 
of the Lorentz invariance |15l EH EH HE which provides also an example of a covariant, power-counting 
renormalizable theory and is represented by a simplest power-law / (R,T, R a pT a P) gravity. 

We have demonstrated that the gravitational held equations in such modified gravity theories admit a 
decoupling property with respect to certain classes of nonholonomic frames, which allows us to generate 
exact solutions for very general off-diagonal forms. The corresponding integral varieties of solutions are 
parameterized by generating and integration functions and various classes of commutative and noncom- 
mutative symmetry parameters. For certain nonholonomic configurations, it is possible to re-dehne the 
generating functions and effective sources of matter fields in such a way that the f(R )-terms are equivalently 
encoded into effective Einstein spaces with complex parametric nonlinear structure for the gravitational 
vacuum. We argue that certain nonholonomic configurations model also covariant gravity theories with 
nice ultraviolet behaviors and seem to be (super-)renormalizable in the sense of Horava-Lifshitz gravity 

[HI [ 20 IHH E21 [13 E3 EH ESI SU SSI 06j 07]. 

Notwithstanding the fact that the various f(R) modified theories and general relativity are actually very 
different theories, the off-diagonal configurations and nonlinear parametric interactions considered in GR 
may encode various classes of such modified gravity effects and explain alternatively observational data for 
accelerating cosmology and certain effects in dark energy and dark matter physics. In both cases, it is 
possible to find cosmological solutions and reconstruct the corresponding action. In the already mentioned 
classes of modified gravity theories with /-modifications [DEEllIMElElElEniEDEElEaElEHlEe} 
E3 E3 El ESI 126] . the dynamics of the matter sector is modeled by a perfect fluid. This is necessary to 
satisfy the continuity equation and guarantee an evolution which is similar to that in GR. For the alternative 
models with nonholonomic configurations [121 E3 EH E21 SU SSI S3 S3 E3 ESI ESI ESI El] , the behavior 
of modified gravity theories is determined by the off-diagonal terms and non-integrable constraints. In 
general, we cannot distinguish the effects of /-modifications from the off-diagonal ones because nonholonomic 
frame transforms mix different classes of nonlinear interactions and parametric constraints. Nevertheless, 
for certain well-defined parameterizations, we can work with effective a(/) and H (/) which, with respect to 
N-adapted frames and for appropriate types of nonholonomic constraints, mimic ACDM cosmology when 
the gravitational background is generically off-diagonal and with a nontrivial gravitational vacuum structure 
(which may be a nonholonomically induced torsion) with an effective cosmological constant. 

Off-diagonal cosmological solutions can be described by a realistic Hubble parameter but with an anoma¬ 
lous behavior for the barionic dark matter as it is shown in the first subsection of Sect. [4] This kind 
of nonlinear parametric evolution allows to reproduce de Sitter like universes modeled on nonholonomic 
backgrounds of certain forms, encoding f(R,T,R a pT a P) gravity. For a corresponding class of generating 
functions, we can model nonholonomic deformations of the ACDM universe with a standard evolution for 
dust matter. It is possible to distinguish corrections with f(T) and/or f{R a pT a terms. The priority of the 
anholonomic frame deformation method (AFDM) is that in such way we can generate analytical and exact 
formulas for the field and cosmological evolution equations, to formulate equivalent modeling criteria, etc. 

Another priority of the AFDM is that we can study the issue of matter instability with various classes 
of modified gravity theories using geometric methods. The equations for the perturbations are complicated 
fourth-order differential equations involving linear perturbations of the Ricci scalar and tensor for different 
classes of linear connections. Nevertheless, we were able to consider specific nonholonomic transforms and 
constraints which allowed us to avoid matter instabilities. In this approach, some viable effective off-diagonal 
Einstein models and f(R) gravities could be elaborated to encode the /(R Qi gT“^) contributions. 

An effective field theory approach to off-diagonal cosmological cosmological configurations can be elabo¬ 
rated in terms of Stiickelberg fields adapted to nonlinear connection structures, see Sect. [5] The constructions 
for generic off-diagonal perturbations of MGTs (in terms of effective speed of sound, gravitational constant, 
Caldwell’s parameter etc.) are linked and confronted with actual observational data. 

We note, finally, that modified gravity theories in general contain ghosts, due to the higher-derivative 
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terms in the action. However, we can select certain ghost-free configurations determined by corresponding 
classes of nonholonomic deformations or constraints. Such models of bi-metric and massive graviton gravities 
were recently studied in mmmmmm- Together with the results in [I5UI6UIT1H8UEU20U2I1 [22] , the 
conclusion is reached that some f(R,T,R a pT a P) models, and their off-diagonal nonholonomic equivalents, 
may possess nice ultraviolet properties and that interesting connections can be established with viable theories 
of quantum gravity. 

Acknowledgments: This work has been partially supported by the Program IDEI, PN-II-ID-PCE- 

2011-3-0256, by an associated visiting research position at CERN, by a DAAD fellowship for Munich and 
Hannover, by MINECO (Spain), grant PR2011-0128 and project FIS2010-15640, by the CPAN Consolider In- 
genio Project, and by AGAUR (Generalitat de Catalunya), contract 2009SGR-994. We thank S. Capozziello, 
N. Mavromatos, S. D. Odintsov, E. Saridakis, D. Singleton, and P. Stavrinos for important discussions and 
support. 


References 

[1] S. Capozziello and V. Faraoni, Beyond Einstein Gravity; A Survey of Gravitational Theories for Cos¬ 
mology and Astrophysics, Fundamental Theories of Physics, vol. 170, (Springer Netherlands, 2011), 467 
PP 

[2] S. Capozziello and M. De Laurentis, Phys. Rept. 509 (2011) 167 

[3] S. Nojiri and S. D. Odintsov, eConf C0602061 (2006) 06 [Int. J. Geom. Meth, Mod. Phys. 4 (2007) 
115] 

[4] S. Nojiri and S. D. Odintsov, Phys. Rept. 505 (2011) 59 

[5] E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70 (2004) 043539 

[6] E. Elizalde, S. Nojiri, S. D. Odintsov and P. Wang, Phys. Rev. D 71 (2005) 103504 

[7] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, JCAP 0502 (2005) 010 

[8] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov and S. Zerbini, Phys. Rev. D 73 (2006) 084007 

[9] G. Cognola, E. Elizalde, S. D. Odintsov, P. Tretyakov and S. Zerbini, Phys. Rev. D 79 (2009) 044001 

[10] E. Elizalde, R. Myrzakulov, V. V. Obukhov and D. Saez-Gomez, Class. Quant. Grav. 27 (2010) 095007 

[11] I. Brevik, E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 84 (2011) 103508 

[12] T. Crofton, P. G. Ferreira, A. Padilla and C. Skordis, Phys. Rept. 513 (2012) 1 

[13] A. de la Cruz-Dombriz and D. Saez-Gomez, Entropy 14 (2012) 1717 

[14] R. Myrzakulov, L. Sebastiani and S. Zerbini, Int. J. Mod. Phys. D 22 (2013) 1330017 

[15] S. Nojiri and S. D. Odintsov, Phys. Rev. D 81 (2010) 043001 

[16] S. Nojiri and S. D. Odintsov, Phys. Rev. D 83 (2011) 023001 

[17] J. Kluson, S. Nojiri and S. D. Odintsov, Phys. Lett. B701 (2011) 117 

[18] E. Elizalde, S. Nojiri, S. D. Odintsov and D. Saez-Gomez, Eur. Phys. J. C70 (2010) 351 

[19] S. Vacaru, EPL 96 (2011) 50001 


34 


[20] S. Vacaru, Gener. Relat. Grav. 44 (2012) 1015 

[21] S. Vacaru, J. Math. Phys. 54 (2013) 073511 

[22] S. Vacaru, J. Geom. Phys. 60 (2010) 1289 

[23] S. Vacaru, Int. J. Geom. Meth. Mod. Phys. 6 (2009) 873 

[24] T. Harko, F. S. N. Lobo, S. Nojiri and S. D. Odintsov, Phys. Rev. D84 (2011) 024020 

[25] S. D. Odintsov and D. Saez-Gomez, Phys. Lett. B725 (2013) 437 

[26] Z. Haghani, T. Harko, F. S. N. Lobo, H. R. Sepangi and S. Shahidi, Phys. Rev. D88 (2013) 044023 

[27] S. Vacaru, J. Phys.: Conf. Ser. 543 (2013) 012021 

[28] S. Vacaru, Int. J. Geom. Meth. Mod. Phys. 8 (2011) 9 

[29] S. Vacaru, J. Math. Phys. 46 (2005) 042503 

[30] S. Vacaru, JHEP 04 (2001) 009 

[31] S. Vacaru and D. Singleton, Class. Quant. Grav. 19 (2002) 3583 

[32] S. Vacaru, Int. J. Mod. Phys. D 21 (2012) 1250072 

[33] S. Bailakos, N. E. Mavromatos, V. A. Mitsou, and M. Plionis, Astropart. Phys. 36 (2012) 7 

[34] A. P. Kouretsis, M. Stathakopoulos and P. C. Stavrinos, Phys. Rev. D86 (2012) 124025; 

[35] M. Hohmann, Phys. Rev. D87 (2013) 124034 

[36] C. Castro, Int. J. Theor. Phys. 51 (2012) 3318 

[37] V. A. Kostelecky, N. Russell and R. Tso, Phys. Lett. B716 470 

[38] C. Lammerzahl, V. Perlick and W. Hasse, Phys. Rev. D86 (2012) 104042 

[39] C. Pfeifer and M. N. R. Wohlfarth, Phys. Rev. D84 (2011) 044039 

[40] D. Kramer, H. Stephani, E. Herdlt and M. A. H. MacCallum, Exact Solutions of Einstein’s Field 
Equations, 2d edition (Cambridge University Press, 2003) 

[41] J. B. Griffiths and J. Podolsky, Exact Space-Times in Einstein’s General Relativity [Cambridge Mono¬ 
graphs on Mathematical Physics] (Cambridge University Press, 2009) 

[42] J. Klusoh, S. Nojiri and S. D. Odintsov, Phys. Lett. B726 (2013) 918 

[43] S. Vacaru, Eur. Phys. J. C 75 (2015) 176 

[44] P. Stavrinos and S. Vacaru, Class. Quant. Grav. 30 (2013) 055012 

[45] S. Vacaru, Int. J. Geom. Meth. Mod. Phys. 11 (2014) 1450032 

[46] T. Gheorghiu, O. Vacaru and S. Vacaru, Eur. Phys. J. C 74 (2014) 3152 

[47] S. Vacaru, Eur. Phys. J. C 75 (2015) 176 

[48] G. F. R. Ellis, Inhomogeneity effects in Cosmology, arXiv: 1103.2335 


35 



[49] K. Godel, Rev. Mod. Phys. 21 (1949) 447-450 

[50] A. Krasinski, Inholomogeneous Cosmological Models , (Cambridge University Press, 1997) 

[51] S. A. Appleby and E. V. Linder, Phys. Rev. D 87 (2013) 023532 

[52] S. Nojiri, S. D. Odintsov and D. Saez-Gomez, Phys. Lett. B681 (2009) 74 

[53] E. Elizalde and J. Q. Hurtado, Mod. Phys. Lett. A 19 (2004) 29 

[54] F. G. Alvarenga, A. de la Cruz-Dombriz, M. J. S. Houndjo, M. E. Rodrigues and D. Saez-Gomez, Phys. 
Rev. D 87 (2013) 103526 

[55] A. D. Dolgov and M. Kawasaki, Phys. Lett. B 573 (2003) 1 

[56] V. Faraoni, Phys. Rev. D. 74 (2006) 104017 

[57] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68 (2003) 123512 

[58] A. Starobinsky, JETP Lett. 86 (2007) 157 

[59] W. Hu and I. Sawicki, Phys. Rev. D 76 (2007) 064004 

[60] S. Nojiri and S. D. Odintsov, Phys. Rev. D 77 (2008) 026007 

[61] S. Capozziello and S. Tsujikawa, Phys. Rev. D 77 (2008) 107501 

[62] A. Appleby, R. A. Battye, Phys. Lett. B654 (2007) 7 

[63] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, Phys. Rev. D 77 (2008) 
046009 

[64] K. Bamba, C. -Q. Geng and C. -C. Lee, JCAP 1011 (2010) 001 

[65] D. B. Kaplan, Five Lectures on Effective Field Theory, arXiv: nucl-th/0510023 

[66] C. Burgess, Annual Rev. Nucl. Part. Sci. 57 (2007) 329 

[67] S. Weinberg, Phys. Rev. D 77 (2008) 123541 

[68] C. Cheung, P. Creminelli, A. L. Fitzpatrick, J. Kaplan and L. Senators, JHEP 0803 (2008) 014 

[69] L. Senatore and M. Zaldarriaga, JHEP 1204 (2012) 024 

[70] P. Creminelli, G. D’Amico, J. Norena and F. Vernizzi, JCAP 0902 (2009) 018 

[71] M. Park, K. M. Zurek and S. Watson, Phys. Rev. D 81 (2010) 124008 

[72] J. K. Boomheld and E. E. Flanagan, JCAP 2012 (2012) 039 

[73] R. Jimenez, P. Talavery, L. Verde, M. Moresco, Climatt, et al., JCAP 1203 (2012) 014 

[74] G. Gubitosi, F. Piazza, and F. Vernizzi, JCAP 1302 (2013) 032 

[75] J. K. Bloomfield, E. E. Flanagan, M. Park and S. Watson, JCAP 1308 (2013) 010 

[76] J. K. Bloomfield and E. E. Flanagan, JCAP 1210 (2012) 039 

[77] A. Starobinsky, JETP Letters 68 (1998) 757 

[78] R. Caldwell, A. Cooray and A. Melchiorri, Phys. Rev. D 76 (2007) 023507 

[79] P. Hofava, Phys. Rev. D79 (2009) 084008 


36 



